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SECTION  I 


INTRODUCTION  AND  SUMMARY 


The  objective  of  the  Thermal  Energy  Storage  Demonstration  Unit  program  was 
the  demonstration  of  a thermal  energy  storage  (TES)  unit  capable  of 
providing  thermal  energy  to  a Vuilleuraier  (VM)  cooler  by  releasing  thermal 
energy  which  had  been  stored  in  the  form  a latent  heat  of  fusion  of  a salt. 
The  requirements  were  to  deliver  1 kW  of  thermal  power  for  1 hour  at  a 
temperature  between  1200°  and  1400°F  to  the  hot  cylinder  of  a Vuilleumier 
cooler.  In  addition,  the  TES  capability  had  to  be  sufficient  to  account 
for  all  thermal  losses  occurring  due  to  limitations  on  the  size  and  weight 
of  thermal  insulation  of  the  TES  unit. 

The  program  was  to  generate  background  data  for  the  design  of  flight  test 
devices  which  can  operate  in  a power  range  of  1 to  3 kW,  while  having  a 
lifetime  of  better  than  5 years.  The  data  were  to  include  costing  for  the 
design,  fabrication,  testing,  and  delivery  of  such  devices,  projected 
efficiency  in  satellite  application,  their  weight,  volume,  and  the  operat- 


ing 

temperature  as  a function  of  load. 

The 

program 

was  organized  into  seven  major  tasks: 

I. 

Analytical  Studies 

II. 

Design 

III. 

Fabrication 

IV. 

Test  Setup 

V. 

Testing 

VI. 

Conceptual  Flight  Device  Parameter  Development 

VII. 

Reporting 

1 


The  total  technical  effort  of  the  program  extended  over  a period  of  12- 
months,  and  was  concluded  with  the  delivery  of  the  tested  TES  unit  mounted 
on  a hot  cylinder  of  a Vuilleumier  cooler. 

For  successful  mating  of  the  thermal  energy  storage  unit  with  the  Vuil- 
leumier cooler,  the  heat  pipe  concept  was  considered  for  transferring  the 
stored  energy  from  the  storage  material  to  the  point  of  energy  requirement. 
This  appeared  most  valid  in  this  application  in  which  the  amount  of  thermal 
energy  is  large  and  the  heat  transfer  surface  of  the  hot  cylinder  over 
which  the  thermal  energy  had  to  be  applied  was  relatively  small. 

As  the  optimization  studies  indicated,  weight  considerations  and  thermal 
loss  considerations  dictated  a unit  with  a relatively  large  diameter,  which 
was  about  one-third  the  length  of  the  unit.  The  largest  acceptable 
diameter  of  the  unit  was  determined  by  the  restriction  that  the  temperature 
at  the  hot  cylinder  of  the  Vuilleumier  cooler  should  not  vary  by  more  than 
50^F  from  the  beginning  of  the  energy  extraction  process  until  the  comple- 
tion of  the  freezing  of  the  thermal  energy  storage  salt.  This  stipulation 
determined  the  largest  thermal  energy  storage  material  thickness  in  the 
unit  and,  thus,  the  greatest  diameter  for  an  annular  design  of  the  storage 
compartment . 

During  the  heat  pipe  design,  the  elevation  pressure  drop,  which  is  deter- 
mined by  the  height  to  which  the  working  fluid  of  the  heat  pipe  has  to  be 
raised,  proved  to  be  the  major  pressure  drop  component.  The  final  wick 
design  of  the  heat  pipe  had  to  utilize  almost  75  percent  of  its  theoretical 
suction  capacity  to  overcome  the  elevation  pressure  drop.  The  wick  design 
proved  to  be  adequate,  although  great  uncertainty  existed  for  a time  when 
the  heat  pipe  did  not  function  properly  for  reasons  which  later  were  traced 
to  a gas  leak.  The  dynamic  behavior  of  wick  structures  has  not  been  fully 
investigated  and  therefore  introduced  uncertainties  into  the  design  of 
large  diameter  heat  pipe  structures  which  have  to  be  operated  in  a 1-g 
environment . 


2 


The  test  results  seemed  to  indicate  that  the  eutectic  salt,  MgF^-LiF-KF, 
used  in  the  thermal  energy  storage  demonstration  unit,  released  only  about 
50  percent  of  the  fusion  energy  of  the  individual  components  of  the 
eutectic  mixture  at  the  eutectic  temperature.  The  total  stored  energy 
could  therefore  not  be  extracted  within  the  design  temperature  range  of 
1250°  -25°F. 
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SECTION  II 


THERMAL  ENERGY  STORAGE  DEMONSTRATION  UNIT  DESIGN 


2.1.  GENERAL  APPROACH 

The  specific  requirements  for  the  thermal  energy  storage  demonstration 
unit  were: 

Thermal  capacity  to  1 kW-hr  thermal 

Power  to  the  hot  cylinder  1 kW  thermal 

Hot  cylinder  operating  tempeature  1250  +25°F 

The  energy  to  the  unit  was  to  be  supplied  from  an  electric  source  which  was 
to  energize  the  thermal  energy  storage  unit  without  the  use  of  electrical 
heaters,  i.e.,  using  the  container  material  of  the  unit  as  resistance 
elements  for  the  electrical  current. 

2.1.1  HEATING  CONFIGURATION 

Considerable  thought  was  given  to  eliminating  electric  heaters  and  using 
solely  the  structural  material  of  the  thermal  energy  storage  unit  as  the 
electrical  pas.;  fwjp  r.-si.stance  heating.  For  a reasonable  current  level, 
the  electric  pass  has  to  be  long  and,  therefore,  many  structural  elements 
of  the  thermal  energy  storage  unit  would  have  had  to  be  placed  electrically 
in  series  with  each  other.  This  would  have  required  electrical  insulation 
of  structural  components  by  insulating  material,  which  appeared  very 
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difficult  to  achieve  when  at  the  same  time  requiring  high  thermal  con- 
ductivity for  maintaining  a reasonable  temperature  distribution  in  the  TES 
unit.  An  even  more  serious  problem  was  thought  to  be  the  galvanic  effect 
to  which  the  thermal  energy  storage  material  would  be  exposed  and  about 
which  little  or  nothing  was  known  for  the  ternary  eutectic  salt.  Any 
galvanic  separation  of  the  individual  elements  could  cause  a serious  attack 
on  the  container  material.  It  was  therefore  suggested  that  the  elec- 
trically resistive  heating  of  the  container  material  should  first  be 
investigated  in  a smaller  device  specifically  developed  for  such  a test. 
This  test  could  furnish  long  time  operating  data  in  parallel  with  the 
operation  of  the  full  scale  thermal  energy  storage  unit  which  could  be 
heated  with  standard  resistance  heaters. 

The  elimination  of  electric  heaters  appeared  to  be  of  great  advantage,  as 

many  electric  heater  failures  have  been  experienced  with  the  VM  cooler. 

When  the  hot  cylinder  is  heated  directly  by  electric  heaters,  the  heat 

2 2 

transfer  area  is  only  7.241  in.  (46.7  cm  ).  For  a power  input  of  1 kV', 

2 

this  indicates  the  relatively  high  power  density  of  21.4  watt/cm  . 

2.1.2  GENERAL  LAYOUT 

After  electing  separate  electric  heating  for  the  power  input  into  the  unit, 
the  most  appropriate  configuration  for  that  heating  scheme  was  considered. 
The  most  critical  considerations  were: 

a.  Lowest  temperature  drop  across  the  thermal  energy  storage 
material 

b.  Smallest  thermal  losses  from  the  unit 

c.  Lowest  weight  of  structural  material 

d.  Lowest  operating  temperature  of  the  heaters 
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The  design  that  evolved  from  these  considerations  was  a tubular  unit  with 
two  interconnected  annular  compartments  for  the  thermal  energy  storage 
material  and  an  annular  well  for  the  resistance  heaters  as  shown  in 
Figure  1. 

The  highest  temperature  drop  across  the  thermal  energy  storage  material  is 
determined  by  the  greatest  thickness  of  the  TES  material  and  the  heat  flow 
rate  per  unit  area.  The  annular  design  provides  large  surfaces  for  a 
relatively  low  heat  flux  as  well  as  a relatively  small  TES  material 
thickness.  By  placing  the  electric  heaters  into  an  annular  well,  the 
heaters  can  give  off  their  power  from  their  entire  surface,  thereby 
operating  at  a lower  power  density  and  corresponding  lower  operating 
temperature  than  a heater  that  transfers  energy  only  from  one  side  while 
having  its  other  side  insulated.  Furthermore,  by  placing  the  heater  inside 
the  thermal  energy  storage  unit,  the  outer  surface  of  the  TES  unit  is  at 
the  lowest  temperature  during  heating  of  the  unit.  This  minimizes  the 
thermal  losses  from  the  unit. 

For  the  initial  layout  of  the  thermal  energy  storage  unit,  the  effective 

2 

heat  transfer  area  was  about  1700  cm  . During  the  time  when  the  heaters 

supply  energy  for  the  operation  of  the  Vuilleumier  cooler  as  well  as  for 

2 

storage,  the  highest  power  density  is  only  1 watt/cm  . This  compares  very 

2 

favorably  with  the  heater  power  density  of  21.4  watt/cm  when  the  hot 
cylinder  of  the  Vuilleumier  cooler  is  heated  directly  with  electric 
heaters.  Despite  the  fact  that  the  heaters  have  to  transfer  power  to  the 
higher  temperature  of  the  thermal  energy  storage  material,  compared  to  the 
operating  temperature  of  the  hot  cylinder  which  is  about  90°F  lower,  the 
heaters  will  actually  operate  at  a lower  temperature  because  of  the 
considerably  lower  power  density  when  energizing  the  thermal  energy 
storage  unit.  Thus,  the  reliability  of  the  heaters  in  the  thermal  energy 
storage  unit  would  be  higher  than  that  of  the  heaters  that  directly  power 
the  hot  cylinder  of  a Vuilleumier  cooler.  This  might  suggest  that  the 
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urgency  for  eliminating  the  heaters  and  for  turning  to  resistance  heating 
the  container  material  of  the  thermal  energy  storage  unit  might  have  been 
eliminated  by  the  thermal  energy  storage  unit  design. 

2.2  DETAILED  DESIGN  APPROACH 

With  the  establishment  of  the  general  approach  to  the  design  of  the  thermal 
energy  storage  demonstration  unit,  the  detailed  design  of  the  unit  required 
the  selection  of  the  thermal  energy  storage  material  and  the  optimization 
of  the  unit  with  respect  to  size  and  thermal  losses. 

2.2.1  THERMAL  ENERGY  STORAGE  MATERIAL 

The  design  of  the  thermal  energy  storage  demonstration  unit  had  to  consider 
the  selection  of  a thermal  energy  storage  material  that  had: 

a.  High  fusion  energy 

b.  Low  volume,  i.e.,  high  density 

c.  Melting  point  appropriate  for  maintaining  the  desired  operating 

temperature  of  the  hot  cylinder 

d.  Compatibility  with  the  structural  material 

e.  Stability  over  many  melting  cycles 

In  the  Solar  Collector  Thermal  Power  System  program  which  preceded  this 
program,  lithium  fluoride  was  selected  as  the  the  thermal  energy  storage 
material.  Lithium  fluoride  has  a melting  point  of  1558°F  and  a heat  of 
fusion  of  0.259  kWh/kg.  For  the  solar  collector  thermal  power  train, 

lithium  fluoride  proved  to  be  the  appropriate  thermal  energy  storage 
material  as  it  had  the  melting  temperature  that  could  provide  the 
Vuilleumier  cooler  with  the  thermal  energy  at  the  desired  temperature  as 
shown  in  Figure  2.  Since  in  the  present  application  the  thermal  energy  did 
not  have  to  be  transferred  aoros.s  a thermal  radiation  joint,  which  required 
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a temperature  differential  of  220°F  in  the  thermal  train,  a thermal  energy 
storage  material  with  a lower  melting  temperature  appeared  to  be  appro- 
priate. A number  of  the  most  promising  thermal  energy  storage  materials 
are  presented  in  Table  1,  while  Figure  3 presents  pertinent  phase 
diagrams.^  From  Table  1 the  binary  eutectic  67  LiF-33  MgF^  and  the  ternary 
eutectic  LiF-30  MgF2~6  KF  were  selected  for  possible  application  in  the 
thermal  energy  storage  unit.  In  Figures  4 and  5 the  calculated  temperature 
distributions  in  the  thermal  energy  storage  unit  are  shown  for  the  two 
materials.  The  ternary  eutectic  LiF-MgF2~KF  appeared  to  provide  the  more 
appropriate  melting  temperature  for  achieving  the  desired  operating  con- 
ditions at  the  hot  cylinder  of  the  Vuilleumier  cooler.  The  inner  surface 
of  the  hot  cylinder,  which  is  cooled  by  the  flowing  helium  gas  in  the 
cooler,  would  have  a temperature  of  1267. 5°F  at  the  beginning  of  the 
thermal  energy  extraction  cycle,  while  at  the  end  of  the  extraction  cycle 
the  temperature  will  have  dropped  to  1227. 5°F.  This  temperature  range  was 
well  within  the  desired  operating  range  of  1250  +25°F.  The  temperature 
level  as  well  as  the  temperature  range  for  a system  using  the  binary 
eutectic  LiF-MgF^  appeared  not  to  be  better  than  that  of  the  ternary 
eutectic. 

Another  ternary  system  which  appeared  to  have  a melting  point  near  the 
desired  operating  temperature  for  the  storage  material  was  62  LiF-19  NaF- 
19  MgF^.  It  fuses  at  708°C  (1306°F).  However,  this  system  was  recognized 
as  unstable  as  it  can  transition  into  a system  which  has  a melting 
temperature  of  only  630°C  (1166°F)  when  NaF  and  MgF^  form  the  double  salt 
MgF^  NaF.  This  is  in  conformance  with  the  double  melting  point  of  the 
binary  system  consisting  of  NaF  and  MgF^,  as  shown  in  the  phase  diagram  of 
Figure  3- 

The  assumed  thermal  energy  storage  capacity  of  the  ternary  eutectic  mixture 
composed  of  biF-MgF^-KF  was  based  on  the  published  data  of  the  heat  of 


1.  "Phase  Diagrams  for  Ceramists,"  American  Ceramic  .Society  (1964). 
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TABLE  1,  PROPERTIES  OF  APPLICABLE  THERMAL  ENERGY  STORAGE  MATERIALS 


Melting 

Point 

Heat  of 

Fusion* 

Eutectic  Molar  Composition 

°C 

kWh// 

kWh/kg 

CaF^ 

1382 

2520 

0. 105 

MgF^ 

1263 

2305 

0.265 

NaF 

996 

1825 

0.603 

0.216 

56  CaF^  + 44  MgF^ 

985 

1805 

KF 

860 

1580 

0.310 

0. 125 

LiF 

848 

1558 

0.464 

0.259 

67  CaF^  + 33  NaF 

813 

1495 

NaCl 

800 

1472 

0.132 

KCl 

776 

1429 

0.098 

79  CaF^  + 21  LiF 

765 

1409 

67  LiF  + 33  MgF^ 

742 

1368 

(0.537) 

(0.238) 

42.5  LiF  + 57.5  MgF^ 

735 

1355 

64  LiF  + 30  MgF^  + 6 KF 

713 

1315 

(0.465) 

(0.219) 

40  KF  + 60  NaF 

710 

1310 

(0.172) 

LiH 

686 

1267 

0.716 

50  KF  + LiF 

650 

1202 

48  LiF  + 52  NaF 

652 

1206 

(0.261) 

NaF  + NaF-MgF^  + LiF  + MgF^ 

630 

1166 

(0.261) 

NaF  + NaF-MgF^  + LiF  + MgF^ 

708 

1306 

62  ZnF^  + 38  LiF 

622 

1152 

* Values  In  parentheses  are  of  questionable  accuracy  and  subject  to 
modifying  definitions. 
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Figure  2.  Temperature  Distribution  Along  the  Solar  Collector  Thermal  Power  Syst 
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fusion  which  were  given  as  9'1'I  joules/gram  in  Ref.  2 and  as  788  joules/gram 
(0.219  kWh/kg)  in  a communication  by  North  American  Philips  Corporation. 
It  appears  now  that  the  heat  of  fusion  of  the  ternary  salt  of  Ref.  2 was  not 
a measured  value  but  had  been  calculated  with  the  assumption  of  "ideality," 
i.e.,  the  heat  of  fusion  of  the  ternary  eutectic  is  the  sum  of  the  heat  of 
fusion  of  its  constituents  when  melting  as  a single  salt  at  the  eutectic 
temperature.  The  latent  heat  quoted  by  the  second  source  was  an  ideal 
value  based  on  the  difference  in  thermal  energy  content  of  the  salt  at 
temperatures  above  and  below  the  eutectic  temperature  and  extrapolated  to 
the  eutectic  temperature.  The  underlying  assumption  must  have  been  that 
the  eutectic  salt  increases  its  thermal  energy  in  a step  function  at  the 
eutectic  temperature  similar  to  a single  component  salt. 

The  test  results  of  this  program  have  shown  that  the  heat  of  fusion  of  the 
ternary  salt  is  only  about  50  percent  the  heat  of  fusion  calculated  under 
the  assumption  of  ideality.  It  should  have  been  obvious  that  ideality  was 
an  inappropriate  assumption  when  looking  at  the  phase  diagram  of  the  salt 
shown  in  Figure  3-  The  phase  diagram  indicates  clearly  the  formation  of  a 
double  salt,  KF-MgF^,  and  transition  lines.  While  the  instrumentation  and 
general  test  setup  of  this  program  was  not  sophisticated  enough  for 
performing  basic  research  on  thermodynamic  properties  of  binary  and 
ternary  salts,  some  general  behavior  and  thermodynamic  properties  could  be 
demonstrated  quite  clearly.  The  results  which  are  presented  in  Section  IV 
of  this  report  proved  the  selection  of  the  ternary  salt  without  the 
availability  of  proven  and  well  documented  thermodynamic  data  to  be  a 
technical  liability. 

The  ternary  eutectic  LiF-MgF^-KF  i.a  not  available  as  a commercial  salt.  It 
had  to  be  prepared  by  combining  its  constituents  in  a retort.  Based  on  the 
molecular  fraction  of  0.64,  0.30,  and  0.06  for  LiF,  MgF^,  and  KF, 


2.  "Summary  and  Selection  of  Inorganic  Salts  for  Application  to  Thermal 
Energy  Storage,"  by  Alina  Borucka,  Report  ERDA  59,  June  1975. 
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respectively,  the  weight  fractions  of  the  eutectic  are  0.4281,  0.4820,  and 

3 

0.0899.  The  specific  weight  of  the  eutectic  was  given  as  2.1233  gram/cm  . 

Preliminary  heat  transfer  calculations  indicated  that  the  losses  from  the 
thermal  energy  storage  unit  would  be  approximately  200  watts,  if  the  unit 
is  insulated  with  Flexible  Min-K  material.  This  predicted  value  proved  to 
be  correct  for  the  steady  state  condition  at  1300°F,  as  will  be  seen  from 
the  test  data  in  Subsection  4.2.1.  The  thermal  energy  storage  unit 
therefore  had  to  be  sized  for  a total  power  release  of  1.2  kW  for  1 hour  of 
operation.  This  corresponds  to  a total  energy  storage  requirement  of 

1.2  kWh  or  4.32  x 10  joules.  When  the  lower  value  of  the  heat  of  fusion  of 
0.219  kWh/kg  or  0.465  kWh/liter  (based  on  the  liquid  density  at  the  melting 
point)  was  accepted,  the  total  weight  of  the  thermal  energy  storage 
material  was  calculated  to  be  5.480  kg  and  the  required  volume  of  the 
material  2.581  liter  (157.481  in.  ).  Thus,  the  total  material  required  for 
the  eutectic  salt  was 

LiF  2.346  kg 

MgF^  2.641  kg 

KF  0.493  kg 

Total  5.480  kg 

2.2.2  OPTIMIZATION  OF  THERMAL  ENERGY  STORAGE  DEMONSTRATION  UNIT 

The  thermal  energy  storage  unit  has  as  its  purpose  the  storing  of  thermal 
energy  in  the  form  of  fusion  energy  of  a salt.  It  is  expected  that  by  using 
a thermal  energy  storage  system,  the  overall  weight  of  the  power  supply  for 
a Vuilleumier  cooler  will  be  lighter  and  that  the  power  system  will  have  an 
equal  or  better  reliability  than  a power  supply  which  stores  energy  in 
electrical  batteries.  For  achieving  this  goal,  the  overall  structure 
containing  the  thermal  energy  storage  material  had  to  be  made  as  small  as 
possible  so  that,  by  itself,  it  will  be  light  and  will  require  the  least 
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amount  of  insulation  for  containing  the  thermal  energy  with  the  lowest 
thermal  losses.  Thermal  losses  from  the  thermal  energy  storage  unit  are 
reflected  in  additional  solar  panel  requirements  as  well  as  a larger 
thermal  energy  storage  unit. 

The  requirement  of  testing  the  thermal  energy  storage  unit  in  a 1-g 
environment  imposed  limitations  on  the  design  which  would  not  exist  if  the 
unit  could  have  been  designed  for  exclusive  space  application  and  would  not 
have  to  be  tested  in  the  laboratory. 

The  general  lay-out  of  the  thermal  energy  storage  unit  was  presented  in 
Figure  1.  It  was  indicated  that  electric  heaters  would  supply  energy  to 
the  thermal  energy  storage  salt  and  that  thermal  energy  for  the  operation 
of  the  hot  cylinder  of  the  Vuilleumier  cooler  would  be  transferred  from  the 
thermal  energy  storage  material  by  heat  pipe  action.  Direct  transfer  from 
the  salt  to  the  hot  cylinder  could  not  be  considered  as  the  heat  transfer 
area  of  the  hot  cylinder  is  very  small  relative  to  the  required  amount  cV 
thermal  energy  salt.  Direct  energy  transfer  would  therefore  not  only 
result  in  a very  high  temperature  differential,  but  would  alo  encompass  a 
large  temperature  variation  during  the  extraction  of  the  energy  from  the 
salt.  The  aim  was  to  find  an  optimum  configuration  with  respect  to  weight 
and  thermal  losses  for  the  basic  design  of  the  unit. 

2.2.2. 1 Overall  Configuration 

Several  parameters  affect  the  weight  and  configuration  of  the  thermal 

energy  storage  unit.  The  most  important  parameters  appeared  to  be  the 

thickness  of  the  thermal  energy  storage  salt,  t^,  the  heater  diameter, 

and  the  wick  thickness,  t , which  are  indicated  in  Figure  6.  For  maintain- 

w 

ing  the  most  uniform  temperature  of  the  hot  cylinder  and  the  heat  pipe,  it 
is  desirable  to  make  the  thickness,  t^,  of  the  thermal  energy  storage  salt 
as  small  as  possibl'^  and  to  achieve  a large  surface  area  and  a small  depth 
across  which  heat  flux  occurs.  The  effect  of  the  salt  thickness  on  the 
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length,  L„„,, , of  the  thermal  energy  storage  unit  Is  shown  in  Figure  7 for 
two  designs  which  differ  solely  with  respect  to  the  size  of  the  heater 
diameter,  and  the  vapor  flow  area,  A^.  A correlation  of  all  the 
dimensional  parameters  of  the  unit  are  shown  in  Figure  8.  The  thermal 
losses  for  the  various  configurations  are  presented  in  Figures  9 through 
16.  It  is  quite  obvious  that,  for  the  lowest  losses,  the  thermal  energy 
storage  salt  should  have  a thickness  of  about  0.8  inch.  Figures  10,  12, 
14,  and  16  also  indicate  that,  after  a certain  insulation  thickness  has 
been  reached,  additional  insulation  will  only  add  weight  to  the  unit 
without  substantially  further  reducing  thermal  losses. 


The  results  of  the  thermal  loss  calculations  are  summarized  in  Figures  17 
and  18.  They  confirm  that  the  vapor  flow  width  has  a minimal  effect  on  the 
thermal  losses  and  that,  therefore,  a minimal  effort  is  required  for 
optimizing  on  this  parameter  in  the  sizing  of  the  unit. 


Based  on  the  parametric  design  study  of  the  thermal  energy  storage  unit, 

the  heat  loss  from  the  unit  appeared  to  be  between  185  watts  '’or  -.n  outer 

insulation  diameter  of  10  inches  and  124  watts  for  an  outer  insulation 

diameter  of  15  inches.  These  losses  compare  with  the  losses  of  a directly 

electrically  heated  hot  cylinder  of  a Vuilleumier  cooler  of  about  50  watts 

as  shown  in  Figure  19  for  the  same  type  insulation.  The  design  that 

evolved  from  the  optimization  study  is  shown  in  Figure  20.  The  length  of 

the  unit  is  about  three  times  the  diameter.  In  this  design,  the  thermal 

2 

load  at  the  heater  surface  is  only  1 watt/cm  . 


2. 2. 2. 2 Insulation  Design 

Prior  to  optimizing  the  insulation  of  the  TE.S  demonstration  unit  by 
employing  a superinsulat ion  consisting  of  radiation  shields  between  layers 
of  conduction  insulation,  it  was  decided  to  use  less  expensive  insulation 
made  of  two  conduction  insulation  materials.  Flexible  Min-K  and  Cerafelt. 
The  insulation  layers  exposed  to  the  highest  temperature  were  to  bo 
Flexible  Min-K.  For  insulating  the  unit  below  600°F,  Cerafelt  was  chosen. 
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Figure  20.  Optimum  Design  for  a 1-kW-Hr  Thermal  Energy  Storage  Unit 


For  the  final  design  of  the  thermal  energy  storage  unit,  the  thermal  losses 
as  a function  of  insulation  thickness  and  outer  insulation  diameter  are 
shown  in  Figure  21.  The  temperature  distribution  through  the  insulation 
for  three  insulation  thicknesses  are  presented  in  Figure  22.  Based  on 
these  parametric  studies,  the  amount  of  insulation  material  for  each  type 
and  the  total  insulation  weight  were  calculated  as  shown  in  Figure  23. 

2. 2. 2. 3 Heat  Pine  Design 

'With  the  basic  dimensions  of  the  thermal  energy  storage  unit  established, 
the  heat  pipe  parameters  had  to  be  determined.  The  analysis  and  heat  pipe 
computer  program  that  had  been  developed  for  the  design  of  the  15-ft  long 
heat  pipe  could  not  be  used  because  of  the  substantial  differences  between 
the  thermal  energy  storage  unit  and  the  heat  pipe  with  respect  to  design 
and  operation.  The  15-ft  long  heat  pipe  had  a distinct  evaporator  and 
condenser  section  which  were  separated  by  a relatively  long  adiabatic 
section.  In  the  thermal  energy  storage  unit,  the  adiabatic  section  is 
missing.  Since  the  evaporator  section  constitutes  the  heat  pipe,  the  flow 
rate  of  the  vapor  as  well  as  of  the  fluid  in  the  wick  vary  along  the  heat 
pipe.  Furthermore,  the  vapor  flow  area  of  the  TES  unit  is  the  combination 
of  an  annulus  and  a pipe,  thus  having  a completely  different  hydraulic 
d;amet«r  than  the  15-ft  heat  pipe.  The  vapor  flow  in  the  thermal  energy 
st.->r a;"’  unit  has  to  account  for  "blowing,"  while  the  flow  in  the  wick  has 

* - .■■.JO*:  ion  flow. 

•,  - .•  ir-' , and  26,  the  pressure  drops  for  three  different  vapor 

: ; - ir*- 1 1<  " rnr.  are  presented.  In  order  for  the  heat  pipe  to  operate 
-ti'  * r;ly,  the  total  pressure  drop  cannot  be  larger  than  the  capillary 
pr<-"- j-'  O'-  suction  pressure  available  from  the  wick  structure.  In  all 

• hr'-"  !e- .gn::  th"  "Icvation  pressure  drop  is  t he  highest  wfiilc  the  pressure 
d-op  in  t.h"  wiok  IS  the  lowest.  The  elevation  pressure  drop  is  the 
lominatin/’’  pr- : .'.ur"  drop  and  is  entirely  determined  by  the  diameter  of  the 
th'-rmil  energy  r.t.orage  unit.  As  was  indicated  earlier,  the  diameter  of  the 


unit  was  predicated  on  the  lowest  thermal  losses  from  the  unit.  The 
optimum  wick  design  was  based  on  the  calculations  shown  in  Figure  27. 
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A bolting  cloth  with  Mesh  70  was  selected  for  the  wick  structure.  This 
design  would  generate  a suction  pressure  about  37  percent  above  that 
required.  As  indicated  by  the  calculations,  a finer  mesh  bolting  cloth 
would  apparently  provide  a larger  safety  margin.  There  is,  however,  an 
uncertainty  whether  the  higher  safety  margin  can  actually  be  achieved 
because  of  limitations  in  maintaining  a consistent  wick  structure  with 
effective  small  pore  sizes  between  the  layers  during  the  manufacturing 
process. 

2. 2. 2. 4 Total  Weight  of  Optimum  Thermal  Energy  Storage  Unit 

The  total  weight  of  the  thermal  energy  storage  unit  is  comprised  of  the 
weight  of:  (1)  the  container  material.  Inconel  600;  (2)  the  wick  filled 
with  the  sodium  working  fluid;  (3)  the  insulation;  and  (4)  the  thermal 
energy  storage  material. 

Container  (Inconel  600) 

Total  Volume 

Specific  Weight 
Weight 

Wick  with  Sodium 

Total  Volume  18.70  inch^ 

Wick  Material  M = 70  inch  \ r 0.00375  inch 

Specific  Weight  of  Wick 

Material  0.302  Ib/inch' 

Effective  Volume  Density  20.33  percent 

Specific  Weight  of  Sodium  0.0279  Ib/inch^ 


1 


47.48  inch^ 

0.302  Ib/inch^ 
14.34  lb 
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TES  PROGRAM  T3  K- 1 .62x10-3xM- .797  CM? 


PROGRAM  23 A 0-002- 000 

PARAMETRIC  ST'JDY  OF  TES  HEAT  PIPE  22  JULY  1975 


WORKING  FLUID  SODIUM, 

22.991 

GRAM/MJL 

OPERATING  TEMPERATURE 

1250.000 

OEG.F 

POWER 

1200.00 

WATT 

heat  pipe  DIAMETERS 

4.776 

INCH  4.3  76  INCH 

.4  00  INC 

WICK  length 

14  .620 

1 NCH 

PROPERTI  ES 

FLUID 

VAPOR 

latent  heat,joules/g 

.401E  04 

surface  TErS  I 0N,DYN7CM 

.134E  03 

viscosity,g/sec-cm 

.195E-02 

.106E- 03 

SPECI FIC  GRAVI TY,G/SEC3 

.789E  00 

.328E-04 

PARTIAL  PRESSURE, DYN/CM2 

.112E  06 

SONIC  VELOCITY,  CM/S  EC  75656.06 


.VALL  WICK  thickness  .0^0  INCH 

WICK  AREA  3.871  Cf.2,  VAPOR  AREA  11.646  CM2,  AREA  RATIO  3.008 
ELEVAU-JN  pressure  9386.91  OYNA;m?,  VELOCITY  784 .04A  CM/SEC 


MESH 

'.VI R E 

PERMEAIML  1 TY 

SUCTI  0 

N DELP  'W 

DELP  V 

DELP  T 

SAFET 

lN-1, 

1 NCH 

CM2 

DYN/CIA2 

DYN/CM2 

DYN/CM2 

DYN/CM2 

PERC 

20 

.009U 

.143E-03 

3630. 

24. 

836. 

1024  7. 

-64  .57 

30 

.0065 

.107E-03 

554  7. 

33. 

836. 

1025  6. 

-45 .91 

40 

.0065 

.854E-04 

8046. 

42. 

836. 

10264  . 

-21.61 

50 

.0055 

.715E-04 

10265  . 

50. 

836 . 

102  72. 

-.07 

60 

.0045 

.618E-04 

12234. 

57. 

836 . 

10280. 

19.01 

70 

.003  7 

.54  7E-04 

14061. 

65. 

836. 

10288. 

36 .68 

80 

.003  7 

.4  92  E- 04 

169  )4. 

72. 

836. 

10295  . 

64  .30 

90 

.0035 

.448E-04 

19556. 

79. 

836. 

10302. 

89.83 

Fif;uro  21.  Ilr>.-it  Pipe  Do.".ir,n  Calculations  for  Final  Doair’;n  of  a 1 -kW-Hr 
Thermal  Kncrjiy  otora/ne  Unit 
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Effective  Specific  Weight 
of  Wick 

Weight 


0.0836  Ib/inch^ 
1.564  lb 


Thermal  Energy  Storage  Material 

Material  (Molecular) 

Material  (Weight) 

Heat  of  Fusion 

Specific  Weight 
Energy  Requirement 

Volume 

Weight 

Insulation 

Material 

Weight 


65  LiF  - 30  MgF^  - 6 KF 

42.81  LiF  - 48.21  MgF^  - 8.99  KF 
99.43  watt-hr/lb 

0.0766  Ib/inch^ 

1200  watt-hr 

157.48  inch^ 

12.07  lb 

Flexible  Min-K  and  Cerafelt 
7.96  lb  (189  watt  loss) 

8.81  lb  (179  watt  loss) 

10.59  lb  (170  watt  loss) 


Based  on  the  above  tabulated  weight  values  of  the  four  major  components  of 
the  thermal  energy  storage  unit,  the  total  weight  of  the  unit  and  the 
specific  power  weights  are  as  shown  in  Table  2. 


TABLE  2.  TOTAL  WEIGHT  OF  THERMAL  ENERGY  STORAGE  UNIT 


Total  Weight  W^ 
(lb) 

35.93 

36.78 

38.56 


Power  Loss  Pj^ 
(watt) 

189 

179 

170 


Sped  fic 
Weight  W^/E 

( Ib/watt-hr ) 

29.94 

30.65 

32.13 


Sped  fic 
Energy  R/W^ 

(watt-hr/lh) 

33.40 

32.63 

31.11 
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The  weight  distribution  for  the  entire  thermal  energy  storage  unit  is  shown 
in  Table  3. 

TABLE  3.  WEIGHT  DISTRIBUTION  OF  THE  THERMAL  ENERGY  STORAGE  UNIT 

Weight  (Ib)  Percent  of  Total 

IP. 07  32.82 

14.34  38.99 


1.56  4.24 

8.81  23.95 

36.78  100.00 

The  values  in  Table  3 indicate  that  the  greatest  impact  on  the  weight  oP 
the  thermal  energy  storage  unit  could  be  achieved  by  a reduction  :n  the 
container  weight.  The  container  weight  is  primarily  determined  by  the 
thickness  of  the  walls,  which  in  the  optimized  design  was  t = 0.065  inch. 
If  by  further  investigation  it  could  be  established  that  thinner  walls 
would  be  safe  for  containing  the  thermal  energy  storage  material  and  the 
working  fluid  of  the  heat  pipe,  considerable  savings  in  overall  TES  unit 
weight  could  be  achieved. 

The  wick  of  the  heat  pipe  contributes  a surprisingly  low  value  to  the 
overall  weight.  Only  4.24  percent  of  the  total  weight  is  attributed  to  the 
wick,  and  the  wick  structure  itself  adds  only  3-12  percent  to  the  weight. 
The  insulation  material  contrib\ates  a .substantial  amount  to  the  overall 
weight,  despite  the  fact  that  it  has  already  been  optimized  by  using  two 
different  materials  according  to  the  temperature  distribution  in  the 
insulation.  The  effect  of  superinsulation  on  heat  losses  and  weight  is 
discussed  in  the  next  section. 
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2. 2. 2. 5 Alternate  Insulation 


For  practical  reasons  it  was  decided  to  employ  Flexible  Min-K  insulation 
with  a layer  of  Cerafelt  for  testing  the  thermal  energy  storage  demonstra- 
tion unit.  The  weight  evaluation  indicated  that  the  insulation  contributed 
about  24  percent  to  the  total  weight  of  the  thermal  energy  storage  unit. 
It  was  concluded  that  the  weight  of  the  unit  and  the  thermal  losses  could 
be  lowered  by  employing  the  insulating  technique  used  for  achieving  the 
lowest  heat  losses  from  a resistojet. 

The  properties  of  high  temperature  thermal  insulation  materials  as  they 
have  been  reported  in  the  literature  are  presented  in  Table  4.  During  the 
development  of  ion  engines  at  XEOS,  in-house  work  on  the  improvement  of  the 
insulation  of  ionizers  was  carried  out.  This  brought  about  the  results 
plotted  in  Figure  28.  The  performance  of  an  insulation  system,  which  was 
called  Linde  -Super  Insulation  and  was  reported  in  the  literature,  has  a 
characteristi''c*'t;/town  in  Figure  29.  This  insulation  system  is  compared  with 
a system  which  v/as  applied  to  the  ion  engine  and  to  a resistojet  in  a 
slightly  different  and  easier  applicable  form.  The  construction  of  this 
insulation  employs  a nominal  0.020  inch  thick  Fiberfrax  paper  separated  by 
0.001  inch  thick  molybdenum  foil.  For  this  construction  the  Fiberfrax 
paper  had  an  effective  thickness  of  only  0.016  inch.  Fiberfrax  is  an 
alumina-silica  fiber  composition  which  is  processed  to  the  Fiberfrax  paper 
with  a binder. 

Using  the  thermal  conductivity  value  determined  from  the  apparent  thermal 
losses  of  the  resistojet,  the  thermal  losses  from  the  thermal  energy 
storage  unit  were  calculated  as  shown  in  Table  5 and  plotted  in  Figure  30. 
For  comparison,  the  thermal  losses  have  also  been  calculated  for  an 
insulation  with  a thermal  conductivity  of  only  9.7  x 10  ^ B- inch/ft^-hr-°F . 
This  insulation  was  mentioned  in  a paper  as  being  available.  The  weight  of 
the  insulation  as  a function  of  the  total  thermal  losses  from  th"  thermal 
energy  storage  unit  are  presented  in  Figure  31  with  the  weight  of  the 
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TABLE  4.  PROPERTIES  OF  HIGH-TEMPERATURE  THERMAL  INSULATION  MATERIALS 


Maximum 

Service 

Temperature 

Density 

Thermal 
Conductivity 
at  1250°F 

Insulating  Material 

(°F) 

(Ib/ft^) 

(Btu/hr-ft  °F 

Thermoflex  Refractory 

2300 

3.0 

1.4 

X 10"^ 

Fiber  Felt 

24.0 

6.4 

X 10“^ 

Micro-Quartz  Felt 

2000 

3.0 

9.5 

X 10“^ 

6.0 

6.9 

X 

o 

1 

Dyna-Quartz  Felt 

2750 

4.5 

10.0 

7.3 

6.3 

X 10"^ 
X 10"^ 

Min-K  2000  (molded) 

1800 

20.0 

2.9 

X 10"^ 

Copper  flakes  dispersed 
in  fibrous  quartz 

1400 

- 

5 X 

10“^ 

Nickel  flakes  dispersed 
in  fibrous  quartz 

1700 

- 

5 X 

10”^ 

Copper  foil  spaced  by 
fibrous  quartz  paper 

1400 

- 

5 X 

10”’* 

Nickel  foil  spaced  by 
fibrous  quartz  paper 

1700 

- 

5 X 

10““ 

Stainless  steel  spaced 
by  ZrO^ 

2000 

- 

3 X 

10““ 

Titanium  spaced  by 
ZrOp  particles 

2000 

_ 

3 X 

10““ 
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HOT  SIDE  TEMPERATURE  T.  ®C 


f 


i 
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TABLE  5.  THEHMAL  LOSSES  AS  FUNCTION  OF  THERMAL  CONDUCTIVITY  OF  INSULA- 
TION MATERIAL  AND  INSULATION  THICKNESS 


INNLR  (JIAMLTLR 

U length  L 

TEMPERATURE  DIFFERENCE  CELT 

1 NCh 

1 rCH 

DEGREE  F 

4.906 

1^)  . 701 

1171.00 

HF.AT  TRA::ifTR 

LUEmu  ENT,['-| 

rtll/i  T2-llfv-DEG.F  9.71E-03 

GUTF.R  LIAMCTEK 

thickness  TI 

POAER  LOSS 

1 NCH 

1 fCH 

KW 

^.406 

.26  0 

.024 

5.^06 

.5  00 

.013 

6.4  06 

.75  0 

.009 

6.906 

1.000 

.007 

7.406 

1 .25  0 

.006 

7. 'Hi  6 

1.5  00 

.005 

R.4  0L 

1 . 76  0 

.00? 

b . 'K'h 

2.000 

.004 

9. 406 

2.25  0 

.OCh 

9.QQ6 

2 .5  00 

.00- 

13.4  06 

? . !-?0 

.003 

10. ''06 

3.000 

.003 

n.4  J6 

3.25  0 

.003 

11. '106 

3 .5  ()0 

.003 

1/  .4  06 

3 . 7p  0 

. CL  3 

INNER  !.’IA'.,ETLR 

C LL  No  1 H L 

TE'.'PERATHRE  ElFFERE’CE  CELT 

1 ilOL 

1 I\CH 

• DEGREE  F 

u.oor, 

15  . 70 1 

1171 .00 

HFAI  IRANSFEk 

CGFFFIS  1 L.M  ,1  -1 

NLH/FT2-R  !-['EG.F  2.76E-02 

UUTEi:  3IAVETFR 

THICK  LESS  TI 

PUOER  L'JSS 

1 N’C  1 1 

INCH 

5.406 

.25  0 

.C'-b 

5.0J6 

.5u0 

.036 

6.4(.  6 

. h C' 

.026 

6."J6 

1 .000 

.020 

7.406 

1.25  0 

.017 

7. '6)6 

1 .5  00 

.OIj 

b .4  06 

1 . 7^  0 

.01' 

0.'>06 

2.000 

.012 

U.4(l6 

2.25  0 

.01  1 

'>.'>06 

2.5  00 

.011 

1 '.40'. 

2 . h 0 

.'1  n 

10  .'106 

3 .OOu 

.OO'i 

1 1 .40f. 

3.25  0 

.00" 

11. ''O'. 

3.5  00 

.02" 

12.401. 

3.7.0 

.OO'j 
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TABLE  5.  THERMAL  LOSSES  AS  FUNCTION  OF  THERMAL  CONDUCTIVITY  OF  INSULA- 
TION MATERIAL  AND  INSULATION  THICKNESS  (Coni' d) 


INN  (3^  01  Ml.,  LICK 

0 length  L 

(E.V.PEKATULE  (MfEEi(E1LC  [£LT 

ll'Cli 

1 ''10  i 

[lELRCE  F 

A.  906 

13.701 

1171.00 

HEAT  TI(A,NSM:u 

CUEl  EIL  1 E\T,B- 

1 NCH/np-l'F-l'EC.F  3 .OOE-01 

OUTER  DIAV.ETfR 

TtllCKNESS  II 

PU'-'-'ER  U5G 

1 iCH 

1 Nl.H 

KW 

6.A06 

.?3  0 

1.229 

‘3.906 

.3  00 

.633 

6.'i06 

.73  0 

.460 

6 . Q06 

1.000 

.364 

7.406 

1.23  0 

.306 

7.QG6 

1.3  00 

0^  .267 

e.406 

1 . 73  0 

.239 

(<.'506 

2.000 

.219 

9.4  06 

2.230 

.202 

9.906 

2.3  00 

.189 

10.4  06 

2 . 73  0 

.179 

10.906 

3.000 

.170 

11.406 

3.230 

.162 

11.906 

3.3  00 

. 136 

1P.406 

3 . 73  0 

.13  0 

*ST(JP*  0 

EOl T LC U \ 

(Dl  T I'FKE 

'TY  1-9 

1 .000  11 . 701,  1 171 

?.000  ,^';»,1v, 

3.O00  0.0U0Y1, 
l.OnO  '■.dOF,  10.701,  1171., 
0.00(1  .,>0,1?, 

.oi'?:, 

7. coo  r.  . 701,  1171., 

( .coo  .?l,10, 

o.OCO  .0, 
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mal  Losses  as  Affected  by 
Amount  of  Insulation 


insulation  composed  of  Flexible  Min-K.  While  the  losses  are  close  to  200 
watts  for  8 pounds  of  Flexible  Min-K,  the  Fiberfrax-Molybdenum  foil  system 
reduces  the  loss  to  less  than  30  watts  for  the  same  weight. 

The  application  of  the  Fiberfrax-Molybdenum  foil  system  is  much  more 
costly.  As  the  insulation  has  to  be  designed  to  minimize  the  heat  flux 
along  the  metal  foil,  the  seams  of  the  insulation  cannot  be  permitted  to 
run  perpendicular  to  the  surfaces.  Therefore,  the  insulation  has  to  be  cut 
to  different  widths  for  each  layer,  which  has  to  be  wrapped  individually 
over  the  thermal  energy  storage  unit.  The  use  of  this  type  of  insulation 
would  have  contributed  little  to  the  demonstration  of  the  thermal  energy 
storage  unit  concept.  In  the  discussion  on  the  testing  of  the  TSS 
demonstration  unit  it  will  become  clear  that  the  relatively  high  thermal 
loss  helped  to  produce  a consistent  and  repeatable  test  condition  for 
determining  the  operation  of  the  TES  unit. 

2.2.3  FINAL  TEIERMAL  ENERGY  STORAGE  DEMONSTRATION  UNIT  DESIGN 

When  placing  the  order  for  the  fabrication  of  the  components  of  the 
optimized  thermal  energy  storage  demonstration  unit,  it  was  impossible  to 
locate  a vendor  who  was  willing  to  fabricate  the  tubes  to  the  design 
dimensions.  The  vendor  who  had  formed  the  tubes  for  the  15-foot  heat  pipe 
was  not  able  to  manufacture  the  smallest  of  the  five  tubes,  which  has  an 
inner  diameter  of  0.493  inch.  He  stated  that  his  manpower  was  no  longer 
available  to  perform  this  work,  although  he  was  willing  to  undertake  the 
forming  of  the  other  four  tubes.  No  other  vendor  could  be  found  to  even 
respond  to  the  request  for  quote. 

Because  of  the  unexpected  difficulties  in  locating  a vendor  for  manu- 
facturing the  tubes,  it  became  obvious  that  the  thermal  energy  storage 
demonstration  unit  would  have  to  be  designed  around  standard  tubes  or  pipes 
made  of  Inconel  600.  In  Figure  32  the  available  Inconel  600  pipes  are 
shown.  Pipes  were  selected  from  this  stock  sheet  that  permitted  the 
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PIPE  SCHEDULES 


-7U 
60  63 


376  0 


2 062 

48J  Y 


S<H.  40S  373  WAU 

34  WT  rti  rr. 


I0*>  103  300  WAil 

S4  T4  ¥<3  fft  ft 


I2*-I<h  105  • 500  WAll 

45  42  W3  H»  M 


ITQ  2 1 

\6'J  1 

1 437 
223  5 i 

1.593 

?4’.  } 

Figure  32.  Available  Inconel  600  Pipes 


manufacture  of  the  thermal  energy  storage  demonstration  unit  with 
some  small  variations  in  the  optimum  design.  The  most  appropriate  pipe 
sizes  for  satisfying  the  design  requirements  were: 


3/8"  nominal 
2"  nominal 
2-1/2"  nominal 
4"  nominal 


0.674"  O.D.  X 0.091"  wall 
2.375"  O.D.  X 0.154"  wall 
2.875"  O.D.  X 0.302"  wall 
4.500"  O.D.  X 0.237"  wall 


Since  a replacement  pipe  for  the  largest  of  the  five  tubes  was  not  available, 
it  still  had  to  be  fabricated  from  available  sheet  metal. 


Figure  33  shows  the  modified  design  that  evolved  from  incorporating  standard 
pipes.  This  design  can  be  compared  with  the  original  design  which  was 
shown  in  Figure  20.  Because  the  smallest  pipe  had  a larger  internal 
diameter  than  in  the  original  design,  the  vapor  flow  passage  in  the 
center  of  the  TES  unit  was  larger  than  required.  This  displaced  some 
of  the  TES  material  which  had  to  be  compensated  for  by  lengthening  the 
unit  approximately  2 inches.  The  lengthening  of  the  unit  resulted  in 
a 3 percent  increase  in  total  weight  as  given  in  the  following  weight 
distribution : 


Design 

Weight,  lb 

Percent 

Original 

TES  material 

12.07 

32.82 

Container  material 
Wick  - screen  material 
- sodium 

1.15  lb 
0.42  lb 

14.34 

38.99 

Insulation 

8.81 

23-95 

Tota  1 

36.79 

100.00 
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56721  A 


Modified 


TES  material 
Container  material 
Wick  - screen  material 
- sodium 
Insulation 

Total 


12.07  31.89 

15.17  110.08 

1.32  lb 
0.48  lb 

8.84  23.28 

37.88  100.00 


The  heater  well  in  the  modified  design  had  a slightly  larger  width.  The 
original  design  called  for  a gap  of  only  0.128  inch.  The  modified  design 
resulted  in  a gap  of  0.136  inch.  With  the  maximum  heater  diameter  of  0.130 
inch,  this  design  was  expected  to  operate  in  the  radiation  heat  transfer 
mode.  Heat  transfer  from  the  heater  to  the  walls  of  the  well  solely  by 
radiation  heat  transfer  had  been  assumed  in  the  original  heat  transfer 
calculations  and  in  the  selection  and  sizing  of  the  heaters.  Any  heat 
transfer  by  conduction  had  been  considered  a bonus  in  the  operation  of  the 
heater,  as  it  would  slightly  decrease  the  heater  operating  temperature. 

When  compared  with  the  original  design,  the  modified  design  resulted  in  a 3 
percent  increase  in  weight.  The  heater  surface  was  slightly  larger  because 
of  the  greater  length  of  the  unit.  The  design  modification  required  for 
using  standard  - though  still  not  easily  available  -Inconel  600  tubes  for 
the  fabrication  of  the  TES  unit  appeared  to  be  acceptable. 

2.2.4  HOT  CYLINDER  DESIGN 


Prior  to  machining  of  the  hot  cylinder,  the  design  was  evaluated  for 
maximum  stresses  and  their  location.  This  was  deemed  important  as  the  hot 
cylinder  was  to  become  an  integral  part  of  the  thermal  energy  storage  unit. 
The  pertinent  dimensions  and  tolerances  of  the  hot  cylinder,  as  shown  in 
Hughes  Aircraft  Company  Drawing  X3180095,  are  repeated  in  Figure  34.  When 
using  the  relations  given  in  Roark:  Formulas  for  Stress  and  Strain,  an1 
the  stress  concentration  factors  given  by  Peterson:  Stress  Concentration 
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Factors,  stresses  as  high  as  157,000  psi  with  the  lower  limit  of  the  plate 
thickness  of  0.070  inch  were  calculated.  The  calculations  took  into 

account  the  curvature  of  the  end  plate  which  attaches  to  the  cylinder  under 
an  angle  of  30.93  degrees  as  shown  in  Figure  3^.  When  using  a simplified 
approach  for  calculating  maximum  stresses  by  matching  of  moments  and 
deflections  of  a curved  plate  with  a cylinder,  stresses  as  high  as  132,000 
psi  were  calculated. 

Data  given  by  the  Carpenter  Technology  Corporation  for  Rene  41 
(Vacumeltrol  4l)  were  correlated  with  the  Larson-Miller  parameter 

= (460  + T)  X (C  + log^Q  t)  X 10“^ 

where : 

T = temperature,  °F 
t = time,  hours 

The  best  fit  was  found  by  making  the  constant  C equal  to  20.  The  cor- 
relation with  the  original  data  is  shown  in  Figure  35.  For  an  operating 
temperature  of  1250°F  and  an  operating  lifetime  of  5000  hours  and  20,000 
hours,  the  Larson-Miller  parameters  are  40.53  and  41.55,  respectively. 
Based  on  the  correlation  of  Figure  35,  the  maximum  allowable  stresses  for 
rupture  are  55,000  and  69,000  psi,  respectively.  For  structural  stability 
and  maximum  allowable  distortion,  the  allowable  stresses  should  be  kept 
below  these  values. 

For  comparison,  the  Larson-Miller  parameter  correlation  with  ‘•he  rapture 
stress  for  Inconel  X750  and  Inconel  718  are  shown  in  Figures  36  and  37. 
The  constant  C is  17.5  for  Inconel  alloy  X750  and  25  for  Inconel  alloy  718. 
The  Larson-Miller  parameters  for  the  two  lifetimes  are  tabulated  below. 
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igure  35.  Rupture  Stress  of  Rone  ')!  (Carpenter  Vaeu"i-3ltrol  Ul) 

(Technical  Data,  Carpenter  Technology  Corporation  6/70) 
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Larson-Miller 

Parameters 

Inconel  X750 

Inconel  718 

/ 

Rene 

5,000 

hours 

36.25 

ii9.08 

no. 53 

20 , 000 

hours 

37.28 

50.11 

41.55 

The  maximum 

allowable 

stresses  for  the 

various  material 

3 are  given 

following  tabulation. 

Maximum  Allowable  Stresses 

Inconel  X750 

Inconel  718 

f 

Rene 

5,000 

hours 

33,000  psi 

50,000  psi 

69,000  psi 

20,000 

hours 

2iJ,000  psi 

25,000  psi 

55,000  psi 

When  the  maximum  allowable  stresses  for  rupture  (this  does  not  consider  the 
maximum  allowable  stresses  for  structural  stability)  were  compared  with 
the  calculated  maximum  stresses  in  the  hot  cylinder  under  the  stated 
operating  conditions,  it  was  felt  that  it  would  be  of  value  to  ascertain 
whether  modifications  could  be  made  in  the  design  without  affecting  the 
operating  characteristic  of  the  hot  cylinder. 

The  most  obvious  improvement  was  the  increase  in  the  plate  thickness  and 
the  tightening  of  its  tolerance.  The  results  of  the  parametric  study  are 
shown  in  Figure  3B.  It  is  interesting  to  note  that  just  the  tolerance  of 
the  plate  thickness  permits  the  maximum  stress  to  vary  between  130,000  and 

96.000  psi.  Without  increasing  the  cylinder  wall  thickness,  which  was 
0.0955  inch  in  the  original  design,  the  plate  thickness  should  be  at  least 
0.09?  inch.  With  this  plate  thickness,  the  maximum  stresses  are  equal  at 

68.000  psi.  A slight  increase  of  the  cylinder  wall  thickness  would  bring 
the  maximum  stresses  down  to  the  allowable  stress  for  20,000  hour  operation 
with  Here  i) l . 
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When  the  manufacture  of  the  hot  cylinder  was  invest If'ated , it  became 
apparent  that  the  availability  of  material  in  the  most  desirable  sizes  was 
very  limited  and  that  the  cost  of  the  material  was  relatively  high.  Since 
the  difficulties  in  obtaining  the  material  increase  with  the  diameter,  it 
became  obvious  that  the  hot  cylinder  should  be  designed  for  the  smallest 
bar  diameter.  The  original  design  called  for  a bar  stock  diameter  of  at 
least  3.875  inches.  This  size  was  needed  for  machining  the  hot  cylinder 
with  the  mounting  flange  as  a single  unit. 

A considerable  saving  in  material  and  machining  could  be  achieved  by 

designing  the  hot  cylinder  with  the  smallest  mounting  flange  possible. 

Calculations  indicated  that,  for  a very  conservative  allowable  bearing 

stress  of  20,000  psi  and  a sheai’  stress  of  15,000  psi,  the  flange  diameter 

had  not  to  be  larger  than  2.292  inches,  and  the  thickness  of  the  flange  not 

more  than  0.029  inch.  A hold-down  ring  with  the  diameter  equal  to  that  of 

the  original  flange,  made  of  plate  material,  could  be  used  for  securing  the 

hot  cylinder  to  the  cooler  body.  Such  design  modification  appeared  to  be 

/ 

especially  useful  if  Rene  4l  was  to  be  used,  as  the  machinabil i ty  of  cnat 
material  is  considered  extremely  poor.  The  proposed  and  accepted  desicn  of 
the  hot  cylinder  is  presented  in  Figure  39  with  the  detailed  design  of  the 
TES  demonstration  unit.  This  change  did  not  require  modifications  to  the 
cooler.  All  design  changes  assured  better  dimensional  stability  and  easier 
fabrication  of  the  hot  cylinder.  Furthermore,  the  material  selection 
became  less  critical  and  made  the  use  of  Inconel  718  possible  without 
reservations  for  the  present  application. 


2.3  SPECIAL  HEAT  TRANSFER  CONSIDERATIONS 


2.3.1  EFFECT  OF  VOLUMETRIC  CHAN'OE  OF  THERMAL  EN'ERilY  .STORAGE  MATERIAL 


With  the  sp.''cific  design  of  the  thermal  ener'gy  .storag-'  iemonr-tration  unit 
e.stabl  ished , the  rr>eondary  ‘■'ffect  associatei)  with  th.’  volumetric  change 
the  thf'rmal  .'tiergy  rd  or- age  matt’riil  during,  melting  and  solidification  w 
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investigated.  In  the  solid  state  the  thermal  energy  storage  material 
occupies  only  76.1  percent  of  the  liquid  volume.  The  container  volume  of 
the  thermal  energy  storage  unit  had  to  be  designed  for  the  largest  volume, 
i.e.,  the  volume  of  the  liquid  material.  During  extraction  of  the  stored 
thermal  energy,  the  salt  contracts  and  voids  form  at  that  time,  so  that  the 
TES  material  will  no  longer  be  in  contact  with  the  entire  container  wall  as 
shown  in  Figures  40,  4l,  and  42.  For  a 24  percent  volumetric  void,  a 20 
percent  surface  detachment  will  occur.  In  the  worst  case  the  temperature 
drop  across  the  thernal  energy  salt  could  increase  by  about  20  percent  for 
a total  of  cT  = 50°F.  Even  with  this  temperature  drop  the  operating 
temperature  of  the  hot  cylinder  would  still  be  above  1200°F. 

2.3.2  COOLir.'G  OF  HOT  CYLINDER  DURING  PERFORMANCE  TESTING 

For  evaluating  the  thermal  energy  storage  demonstration  unit  performance, 
the  hot  cylinder  was  cooled  with  a flow  of  air  simulating  the  heat  extrac- 
tion from  the  hot  cylinder  wall  by  helium  during  operation  of  the  Vuil- 
leumier  cooler.  The  arrangement  for  testing  the  thermal  energy  storage 
demonstration  unit  is  shown  in  Figure  43.  A plug  shown  in  Figure  44  was 
machined  from  stainless  steel  that  provided  a 0.020-inch  flow  pass  along 
the  hot  cylinder  wall.  The  predicted  heat  transfer  as  function  of  flow 
rate  is  shown  in  Figures  45  and  46.  With  the  0.020  inch  flow  passage,  the 
design  power  is  extracted  at  a flow  rate  of  8.2  x 10  ^ Ib/sec  (3-72  g/sec). 
At  this  flow  rate  the  flow  is  well  in  the  laminar  flow  regime,  and  the  heat 
transfer  varies  relatively  little  with  flow  rate,  i.e.,  the  total  heat 
transfer  changes  by  only  5 percent  for  a flow  rate  variation  of  about  24 
percent . 

The  length  of  the  flow  gap  is  L = 1.025  inch.  With  a flow  gap  of  only 
t = 0.020  inch,  the  entrance  regime  is  very  short  compared  to  the  entire 
flow  length.  For  evaluating  the  Nusselt  number,  the  relation  given  by 
Eckert  and  Drake  in  Heat,  and  Mass  Transfer,  2nd  Edition,  p.  107,  is  used. 
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Figure  40.  Configuration  of  Thermal  Energy  Storage  Demonstration  Unit 
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Volumetric  Void  as  a Function 
of  Half  An^le  for  Thermal 
Energy  Storage  Demonstration 
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The  non-dimensional  parameter  l/(Re^  Pr  d/x)  is  about  0.05  so  that  the 
average  effective  Musselt  number  is  about  5 for  the  case  under  considera- 
tion. As  shown  by  the  correlation  given  by  Eckert  and  Drake,  Prandt  and 
various  other  literature  sources,  the  heat  transfer  is  not  entirely 
dependent  on  the  Nusselt  number  but  is  also  a function  of  several  addi- 
tional parameters.  Inlet  conditions,  roughness  factors,  and  surface 
conditions  are  contributing  parameters  which,  however,  could  not  be 
expressed  mathematically  for  inclusion  in  the  heat  transfer  calculation. 

The  final  temperature  of  the  air  that  cools  the  hot  cylinder  of  the 
Vuilleumier  cooler  and  has  an  inlet  temperature  of  70°f,  is  shown  in 
Figure  47.  The  calculated  pressure  drop  along  the  flow  pass  is  shown  in 
Figure  48.  With  the  expected  pressure  drop  being  approximately  1 psi , it 
was  possible  to  arrange  the  flow  metering  system  to  operate  in  the  chocked 
orifice  configuration.  With  the  additional  pressure  drops  in  the  lines  and 
the  fittings  leading  to  and  from  the  flow  insert,  the  total  back  pressure 
was  estimated  to  be  not  higher  than  10  psig.  The  sonic  orifice  was 
therefore  selected  to  operate  with  an  upstream  pressure  of  about  50  psia  or 
35  psig. 

For  spanning  the  entire  range  of  flow  rates  which  might  be  desirable  for 
testing  with  shop  air  supply,  two  sonic  orifices  were  calibrated,  No.  55-1 
and  No.  47.  Orifice  No.  55-1  permitted  lowering  the  flow  rate  to  less  than 
0.7  g/sec  of  air,  while  with  orifice  No.  47,  the  upstream  pressure  was 

_3 

a'jout  69  psig  at  the  design  flow  rate  of  8.2  x 10  Ib/sec.  The  orifices 
were  calibrated  with  a Wet  Test  Meter  (American  Meter  Company  Mod^l  AL  22), 
which  has  a maximum  deviation  of  -1  percent  from  the  calibration,  and  with 
a 1/4  percent  Dourdon  tube  type  pressure  gage  (Glassco). 

In  Figure  46  tv;o  regimes  are  indicated.  Below  1.36  g/se^'  th--  am  f 
energy  extracted  from  the  hot  cylin<lcr  will  be  equal  t.o  the  iu-a*  - ir  i''  * v 
of  the  air  bet.wecn  the  inlet  tcmpt^rature  and  the  operating  ; .j.  ..  f 


69 


Figure  '48.  Pres.sure  Drop  Along  the  Heat  Transfe 
Path  of  the  Hot  Cylinder  of  the 
Vuilleumier  Cooler 


hot  cylinder,  as  the  air  will  be  heated  to  the  operating  temperature  of  the 
hot  cylinder  when  flowing  through  the  flow  gap.  Above  1.36  g/sec  the 
cooling  air  will  not  reach  the  operating  temperature  of  the  hot  cylirider. 
For  this  regime  the  power  extraction  is  determined  by  the  heat  transfer  to 
the  cooling  air  under  laminar  flow  conditions. 

A testing  problem  was  anticipated  for  the  case  when  trying  to  raise  the 
extraction  rate  to  about  twice  the  design  rate.  As  can  be  seen  fro.m 
Figure  U5,  a heat  transfer  rate  corresponding  to  1?0  to  200  percent  of  the 
design  extraction  rate  lies  in  a regime  in  which  the  flow  is  unstable  for 
the  0.020-inch  gap.  Although  the  0.020-inci  gap  is  equal  to  the  gap  of  t'ne 
hot  cylinder  when  mounted  on  the  Vuilleumier  refrigerator , it  was  antici- 
pated to  increase  the  gap  by  remachining  the  flow  insert  at  a later  date. 
This  would  have  permitted  testing  at  heat  transfer  rate."^  above  the  design 
rates  for  the  device.  The  correlations  in  Figure  '19  indicate  that,  ■.•;ith  a 
flow  gap  of  0.100  inch  the  hot  cylinder  can  be  tested  from  bO  percent  of  the 
design  heat  transfer  without  change  in  flow  regime;  i.e.,  the  entire  flow 
regime  is  turbulent.  Very  consistent  test  data  should  have  been  possible 
with  this  gap  configuration.  The  pressure  drops  that  were  anticipated  for 
that  conf iguration  are  shown  in  Figure  50. 
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SECTION  III 


FABRICATION  OF  THERMAL  ENERGY  STORAGE  DEMONSTRATION  UNIT 


3. 1 MACHINING  AND  INITIAL  ASSEMBLY 

The  thermal  energy  storage  unit  was  fabricated  from  Inconel  600  sheet  metal 
and  Inconel  600  pipes  according  to  the  availability  of  material  and  the 
capability  of  manufacturers  to  form  the  required  components.  Figure  51 
shows  the  machined  parts  in  an  exploded  view  of  the  therm.al  energy  storage 
demonstration  unit  with  the  hot  cylinder  and  the  hold  down  flange.  The 
four  smaller  of  the  five  tubes  were  machined  from  pipes,  while  the  largest 
was  roled  from  0.065-inch  Inconel  600  sheet  metal.  The  same  material  was 
used  for  machining  the  seven  flanges  shown. 

As  indicated  in  a prior  section  of  this  report,  the  original  TES  design  had 
to  be  modified  to  permit  m.anufacture  of  the  tubes  from  available  Inconel 
600  pipes  because  the  tube  forming  vendor  did  not  want  to  bid  on  forming 
the  smallest  of  the  five  required  tubes. 

The  cost  of  manufacturing  the  four  tubes  from  standard  pipes  instead  of 
from  sheet  metal  was  not  lower,  as  one  might  have  expected.  Material  cost 
was  considerably  higher  because  only  thick-walled  pipes.  Schedule  UO,  are 
being  manufactured  and  stocked  by  the  suppliers.  The  cost  of  the  pipt.a  is 
based  primarily  on  the  material  weight  of  Inconel  600,  which  is  relatively 
expens i v" . 

The  hot  cylinder  of  the  Vuilleumier  cooler  was  machined  from  Inconel  7l8, 
2l-inch  diamet.er  tiar  stock.  The  inside  of  Uie  hot  cylindi'r  was  fii’st  rougii 
machine'.  r"moving  tiie  m.ajor  amount  of  mat-Tial.  The  final  dimensions  were 
obtained  by  the  Eloy  process.  The  split  hold-down  flange  was  fabricated 

frem  iinl-'.-,;;  steel. 
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Figure  51.  Exploded  View  of  Thermal  Energy  Storage  Demonstration  Unit 


The  assembly  of  the  thermal  energy  storage  unit  was  initiated  by  electron 
beam  welding  all  joints  required  for  forming  ttie  container  of  the  thermal 
energy  storage  material.  Since  the  heater  well  walls  are  part  of  the 
thermal  energy  storage  material  container,  the  placement  of  the  six 
electric  heaters  had  to  be  included  in  the  initial  asrambly  of  the  unit. 
The  six  electric  heaters  were  first  wound  tightly  over  a 2-inch  diameter 
tube  and  then  transferred  as  a unit  over  the  inner  surface  of  the  heater 
well.  The  outer  tube  of  the  heater  well  was  then  slipped  over  the  heaters. 
The  installation  of  the  heaters  proved  to  be  more  difficult  than  had  been 
anticipated.  The  1 /8-inch  diameter  heaters  could  not  be  wound  uniformly 
enough  to  be  easily  accommodated  even  in  the  0.132  inch  nominal  wide  well. 

With  the  thermal  energy  storage  material  container  mechanically  assembled 
anl  the  heaters  in  place,  the  remaining  welds  were  completed.  At  this 
stage  of  the  assembly,  the  joints  were  checked  with  a helium  leak  detector. 
One  joint  was  found  faulty  and  had  to  be  redone.  This  joint  had  to  be 
welded  with  the  heaters  in  place  in  the  heater  well.  The  extensions  of  the 
six  heaters  caused  the  electron  beam  to  be  deflected  and  to  wander  away 
from  the  joint.  During  weld  repetition,  the  heater  ends  were  .■shielded, 
thus  preventing  the  interaction  of  the  heater  ends  with  the  electron  beam. 

The  thermal  energy  unit  as  it  appeared  at  this  stage  of  the  manufacturing 
procerus  and  prior  to  its  filling  with  the  thermal  energy  storage  material 
is  shown  in  Figure  52. 

3.2  LOADT’IO  OF  THFRMAl,  FWF.RCY  STORAGE  UNIT  WITH  TFS  MATFRIAL 
3.2.)  EVAi.UATIOW  OF  LOADING  PHOCFDURES 


During,  the  initial  de.slgn  of  the  thermal  cnerg.y  steragv;  demonstration  iin  . t . 
it  war  planned  to  fill  t.ho  unit  with  the  thermal  energy  storag.e  rait  by 
first  card,inr,  it  in  a car, ting  form  male  of  Mullite  t.ut'ing.  The  ".ull.te 
tubin'’;  war  acquired  earlv  In  tie'  progr- im  with  Mie  seiecte.i  tube  -.e;- 


reflecting  the  optimum  thermal  energy  storage  unit  design.  When  the 
initial  thermal  energy  storage  unit  design  had  to  be  modified  to  reflect 
the  availability  of  material,  no  standard  Mullite  tubing  could  be  found 
that  would  permit  the  casting  of  the  thermal  energy  storage  salt  to  the 
desired  sizes. 

When  several  alternatives  to  filling  the  unit  with  its  thermal  energy 
material  were  considered,  the  melting  of  the  ternary  salt  directly  in  the 
unit  was  investigated.  The  correct  amount  of  the  three  reagent  grade, 
chemically  pure  salts,  LiF,  MgF.,,  and  Kl'' , which  are  available  in  powder 
form  would  be  poured  into  the  unit  in  their  correct  combination  and  the  six 
heaters  would  be  used  for  melting  the  mixture  directly  in  the  unit. 

To  evaluate  the  casting  process,  a small  container  was  made  of  Inconel  600 
tubing.  The  volume  of  the  container  was  sized  to  hold  100  grams  of  the 
eutectic  mixture  in  powder  form  which  occupies  about  three  times  the  volume 
of  the  eutectic  in  solid  form.  One  hundred  grams  of  the  salt  mixture  in 
powder  form  was  measured  to  occupy  92  milliliters. 

The  container  was  surrounded  by  a shell  heater  which  was  controlled  with  a 
variac.  Two  thermocouples  were  placed  at  the  bottom  and  along  the  side  of 
the  container.  The  output  of  the  thermocouples  was  read  on  a digital 
voltmeter  and  also  continuously  displayed  on  a chart  recorder.  The  results 
of  a typical  melt  are  shown  in  Figure  53.  indicating  the  smallest  rate  of 
increase  in  temperature  at  1310°F  during  the  melting  process.  After  the 
rate  of  increase  in  the  temperature  had  increased  again  for  6 minutes  the 
heaters  were  shut  off. 

After  molting  the  mixture  of  the  three  components  of  the  eutectic,  LiF, 
MgF^,  and  KF,  the  occupied  volume  of  the  container  was  only  one-third  of 
its  total  volume.  Successive  additions  of  the  powder  mixture  to  the 
initial  melt  should  have  made  it  possible  to  fill  the  eont.ainer  with 
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Figure  53-  Temperature  of  Thermal  Knergy  Storage  Material  Capsule  During  Melt-down 


300  grams  of  eutectic  salt.  This  proved  to  be  correct.  However,  dur 
the  filling  process,  the  surfaces  of  the  container  at  the  upper  edges  w 
considerably  contaminated. 

Inspection  of  those  surfaces  indicated  that  if  welding  would  have  t 
intended  without  danger  of  inclusion  of  contaminants  remachining  wc 
have  been  required.  The  test  indicated  the  need  for  rework  of  the  TES  u 
after  filling  with  the  TES  material,  if  the  unit  would  have  been  filled 
casting  the  material  in  the  unit  itself. 

As  the  risk  of  damaging  the  hardware  by  the  above  outlined  filling  apprc 
appeared  to  be  high,  it  was  decided  to  perform  the  casting  in  a cast 
fixture  as  had  originally  been  intended. 

3.2.2  CASTING  OF  TES  MATERIAL 

For  the  casting  of  the  thermal  energy  storage  salt,  a casting  fixture 
machined  from  graphite  (grade  AGSX)  which  was  secured  from  Union  Cart 
Corporation.  The  casting  fixture  was  sized  for  casting  one-third  of 
total  thermal  energy  storage  salt  in  each  heat  in  a brazing  retort, 
casting  fixture  design  is  shown  in  Figures  54  and  55,  while  the  fi 
machined  parts  of  the  fixture  are  shown  in  Figure  56. 

The  retort  could  be  heated  with  an  inert  gas  or  in  a hydrogen  atmosphe 
An  argon  atmosphere  was  selected  for  the  casting  to  prevent  any  chemi 
reaction  with  the  three  salts,  LiF,  Mg^F,  and  KF.  The  temperature 
raised  to  1350*^F,  which  is  abo\it  40°F  above  the  melting  temperature 
ternary  eutectic,  and  was  held  at  this  temperature  for  about  1 hour 
ensure  full  melting. 

A total  of  four  casts  were  made  to  account  for  a possible  failure  of  a c 
and  to  assure  suffici<^nt  amount,  of  salt  available  for  filling  the  TES  i. 
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after  the  casting.  This  proved  to  be  a very  prudent  precaution  as 
cast  turned  out  to  be  short  by  sevei’al  grams.  P’urthcrmoro , the  res' 
pieces  of  the  salt  showed  voids  in  the  center.  The  salt  appeared  tc 
solidified  from  the  bottom  as  well  as  from  the  top,  leaving  the  c 
portion  full  of  voids  between  a large  solid  bottom  piece  and  a sn 
solid  top  piece.  This  could  have  been  avoided  by  forcing  cooli 
proceed  from  the  bottom  of  the  cast,  which  may  iiave  been  achiev 
forcing  the  inert  gas  around  the  bottom  of  the  test  fixture  durin 
cool-down  and  by  partially  insulating  the  top  part  by  increasing  the 
of  the  graphite. 

The  casting  process  of  the  thermal  energy  storage  material  in  this  pr 
varied  from  the  process  used  by  General  Electric  in  an  earlier  progra 
casting  single  solid  pieces.  At  General  Electric,  the  m.aterial  was  ca 
a vacuum  with  shell  heaters  surrounding  the  single  tube.  After  the 
had  become  molten,  the  casting  fixture  was  lowered  with  the  he 
remaining  in  place.  Thus,  the  lower  portion  of  the  fixture  was  co 
while  the  upper  portion  remained  hot  and  no  voids  could  form.  Thi.a 
relatively  expensive  casting  process  as  a vacuum  chambe:’  was  contami 
by  the  very  corrosive  vapors  of  the  salt  and  had  to  be  thoroughly  cl 
after  its  use. 

3 . 3 FINAL  A.SSEMRL,Y  OF  THE  THERMAL  ENERGY  STORAGE  DEMO:JSTRATION  Ul.'IT 

After  the  thermal  energy  storage  demonstration  iinit  had  been  loaded 
the  cast  TK.G  matr>rial,  the  final  two  welds  of  th"  container  were  compl 
Unfortunately,  at  this  .stage  of  the  manufacturing  proc'uss  it  wa: 
possible  to  leak  check  the  final  t.wo  wehis.  F^d  lowing  the  iw-  w 
wicking,  material,  con.’ i rt i ng,  ol'  s''veral  layerr.  of  bolting  cloth  of  I’O 
was  applied  to  tiie  inner  surface  of  the  outa'r  tube  of  TES  unit 

outer’  .surface  of  t.h'*  storage  enni.  nn<>r,  t Iv’  f 1 ang/'s- , an  1 * h*’-  hot  cyli 
Tire  d i st  r’ iliut  i on  of  th"  wick  mat'-riil  is  .shown  in  Figure  ty.  Whcr’’V’ 
wicking  mat«'rial  was.  covering,  a iieat  t.ransfer  surface,  the  thickness  c 
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wick  was  made  not  larger  than  four  layers  of  bolting  cloth.  All  other 
surfaces  were  covered  with  eight  layers  of  bolting  cloth.  The  appearance 
of  the  thermal  energy  storage  unit  with  all  internal  surfaces  covered  with 
the  wicking  material  is  shown  in  Figure  58.  The  locations  of  the  eight 
pads  which  bridge  the  outside  wick  with  the  wick  of  the  TES  container  are 
evident  in  this  picture. 

The  assembly  of  the  thermal  energy  storage  unit  was  completed  with  the 
welding  of  the  outer  tube  to  the  hot  storage  container.  The  completed 
thermal  energy  storage  demonstrated  unit  is  shown  in  Figure  59. 

Sixteen  nickel-chromium  versus  nickel-aluminum  thermocouples  were  placed 
on  the  inside  and  the  outside  of  the  thermal  energy  storage  unit.  The 
locations  of  these  thermocouples  are  shown  in  Figure  60.  The  outside 
thermocouples  were  attached  to  the  surfaces  by  spot  welding  stainless  steel 
ribbons  to  the  Inconel  over  the  thermocouple  naked  junctions.  The  thermal 
energy  unit  was  then  wrapped  with  five  layers  of  3/8-inch  Flexible  Min-K 
insulation  and  one  layer  of  1/2-inch  8 pound  weight  Cerafelt  insulation. 
The  insulation  was  completed  with  one  wrap  of  aluminum  foil.  The  predicted 
thermal  losses  for  the  thermal  energy  storage  demonstration  unit  with  this 
type  of  insulation  are  shown  in  Figure  6l.  The  completely  insulated  TES 
unit  ready  for  filling  of  its  heat  pipe  with  sodium  is  shown  in  Figure  62. 

3.14  FILLING  OF  HEAT  PIPE  WITH  WORKING  FLUID 

Based  on  the  modified  thermal  energy  storage  demonstration  unit  design,  the 

3 

theoretical  total  wick  volume  was  calculated  to  be  27.383  inches  . For  a 
wick  constructed  of  a 70  Mesh  wire  screen,  woven  with  a 0.00375  inch 
diameter  wire,  the  void  volume  that  is  filled  with  sodium  comprises  79.67 

3 

percent  of  the  total  wick  volume.  Thus,  a sodium  volume  of  21.82  inches 
was  indicated  for  filling  the  heat  pipe  section  of  the  thermal  energy 
storage  unit.  At  the  operating  temperature  of  1250°F,  the  specific  weight 

3 

of  sodium  is  O.7B96  g/cm  , which  determined  the  total  weight  of  the  fill  to 
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INSULATION; 

5 LAYERS  FLEX  MIN-K  3/8-INCH  \ 
C^' THERMOCOUPLE  LOCATION  (BOTTOM)  , ^^yER  CERAFELT  1/2-INCH  8LB.  I 

O THERMOCOUPLE  LOCATION  (TOP) 


V THERMOCOUPLE  LOCATION  (SIDE) 


Figure  60.  Insulation  and  Thermocouple  Locations  of  Thermal  Energy  Storage 
Demonatrat Lon  Unit  for  Vuilleumier  Cooler  Application 
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be  282.3  grams.  Since  the  sodium  was  to  be  filled  at  a temperature  of 

320°F,  at  which  the  specific  weight  of  the  sodium  is  0.913  gm/cm^,  a 

3 

calibrated  fill  volume  of  18.523  inches  was  needed. 

A fill  volume  was  built  from  a stainless  steel  tube  which  had  an  O.D.  of 
1.250  inch,  a wall  thickness  of  0.065  inch,  and  a length  of  I8.8OI  inch. 
The  tube  was  weighed  prior  to  and  after  filling  it  with  sodium.  The  total 
amount  of  sodium  in  the  tube  was  slightly  more  than  282.3  grams  (i.e., 
284.8  grams)  because  of  the  additional  volume  in  the  1/4-inch  tubes  at  the 
top  and  the  bottom  of  the  fill  volume.  After  transferring  the  sodium  into 
the  thermal  energy  storage  unit,  the  weight  of  the  fill  volume  indicated 
that  the  heat  pipe  section  had  been  filled  with  281.1  grams  of  sodium. 

It  was  attempted  to  fill  the  heat  pipe  section  fairly  accurately  with  the 

calculated  required  amount  of  sodium,  although  realistically  the  exact 

amount  of  sodium  needed  could  not  be  verified  until  the  thermal  energy 

storage  unit  had  been  operated.  Of  all  the  individual  wick  volumes  which 

comprise  the  total  volume  of  the  sodium,  the  eight  layers  of  screen  on  the 

inner  surface  of  the  outer  tube  is  the  largest.  Its  calculated  volume  is 

3 

15.631  inches  comprising  57.1  percent  of  the  total  calculated  wick  volume 
alone.  As  it  is  difficult  to  achieve  a perfect  wick  of  consistent  void 
volume  with  bolting  cloth,  the  required  amount  of  sodium  for  filling  the 
the  wick  without  surplus  or  dryness  in  certain  areas  has  to  be  established 
by  testing  the  heat  pipe. 

With  the  filling  of  the  heat  pipe  with  the  working  fluid  the  fabrication 
tasks  were  completed  and  the  thermal  energy  storage  unit  was  ready  for 
testing. 
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SECTION  IV 

TESTING  OF  THERMAL  ENERGY  STORGE  DEMONSTRATION  UNIT 


4.1  TEST  SETUP 

VJith  completion  of  its  fabrication  the  Thermal  Energy  Storage  Demonstra- 
tion Unit  was  integrated  into  the  test  setup  which  is  shown  in  Figure  63- 
The  thermal  energy  storage  demonstration  unit  had  been  supplied  with  six 
individual  heaters  grouped  in  sets  of  two.  The  two  heaters  in  each  set 
were  operated  in  parallel  with  each  other,  with  their  power  being  control- 
led by  a variac.  The  average  resistance  of  each  heater  was  18  ohms,  so  that 
at  an  applied  voltage  of  110  volts  the  power  input  was  about  672  watts  per 
heater  for  a possible  total  power  input  of  about  4000  watts  with  all  six 
heaters.  Thus,  even  with  only  two  heaters  operating,  it  was  possible  to 
charge  the  thermal  energy  storage  unit  within  one  hour. 

The  three  variacs  could  be  adjusted  and  their  power  output  determined  by 
three  sets  of  1 percent  full  scale  accurate  ammeters  and  voltmeters.  For 
controlling  the  temperature  of  the  thermal  energy  storage  unit,  a tempera- 
ture controller  was  used  to  turn  off  the  power  to  the  unit  if  the  preset 
temperature  was  exceeded  and  to  turn  on  the  power  when  the  temperature  of 
the  unit  dropped  below  a preselected  temperature. 

While  the  temperature  of  all  sixteen  thermocouples  could  be  determined  from 
the  output  of  a digital  voltmeter,  twelve  of  the  temperatures  were  recorded 
continuously  on  a 24  point  chart  recorder.  The  gas  flow  for  cooling  the 
hot  cylinder  was  controlled  by  the  upstream  pressure  of  a sonic  orifice, 
and  by  observing  the  pressure  downstream  of  the  orifice  to  assure  that,  at 
all  times  during  a tost,  sonic  conditions  at  the  orifice  would  exist.  The 
pressures  were  measured  with  1/2  percent  accurate  Bourdon  pressure  gages 
which  can  be  seen  in  Figure  64. 
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Figure  63 


Thermal  Energy  Storage  Demonstration  Unit  Test  Setup 
(View  I) 


Figure  6»l,  Thermal  Energy  Storage  Demonstration  Unit  Test  Setup 
(View  II) 


92 


n.2  INITIAL  TESTING 

After  completion  of  the  test  setup  and  the  instrumentation  of  the  thermal 
energy  storage  unit,  testing  of  the  unit  was  initiated.  The  first  objec- 
tive was  to  measure  the  steady  state  thermal  losses  from  the  unit  as  a 
function  of  the  operating  temperature,  and  to  determine  the  effective 
thermal  capacity  of  the  unit,  (W 

4.2.1  STEADY  STATE  THERMAL  LOSSES  AND  THERMAL  CAPACITY 

4. 2. 1.1  Time  Constant  and  Steady  State  Thermal  Losses 

The  thermal  energy  storage  unit  comprises  the  main  structure,  the  wick 
material,  the  heat  pipe  working  fluid,  sodium,  and  the  thermal  energy 
storage  salt.  The  effective  thermal  capacities  of  the  four  components  are: 


Structural  Material  (Inconel) 

1 .865 

E/°F 

27.80? 

Thermal  Energy  Storage  Salt 

4.459 

B/°F 

66.48? 

Sodium 

0. 190 

B/°F 

2.88? 

Wick  Material 

0.193 

B/°F 

2.83? 

Total 

6.707 

B/°F 

100.00? 

7.076 

X 10^ 

joules/°F 

During  the  solidification  of  the  thermal  energy  storage  salt,  a total 
amount  of  thermal  energy  of  4.32  x 10^  joules  was  expected  to  be  released. 
It  was  predicted  that  the  temperature  of  structural  components  of  the 
thermal  energy  storage  unit  would  decrease  from  1275°  to  1235°F,  while  the 
average  temperature  of  the  thermal  energy  storage  salt  would  decrease  from 
1310°  to  1270°F.  The  contribution  of  the  sensible  heat  of  the  thermal 
energy  storage  unit  would  therefore  be 

E ...  = 4.459  X 40  + 2.248  x 50  = 290. 76B 

sensible  ' , 

= 0.307  X 10  joules 
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This  compares  with  the  assumed  total  latent  heat  of  the  tfiermal  energy 
storage  salt  of 

E = M.32  X 10^  joules, 

latent 

i 

The  contribution  of  the  sensible  heat  would  thus  be  only  6.63  percent  of  j 

the  total  thermal  energy  released.  i 

I 

The  steady  state  heat  losses  from  the  thermal  energy  storage  unit  were  to  i 

be  obtained  by  correlating  the  power  input  with  the  steady  state  tempera-  j 

f 

ture  of  the  unit.  For  predicting  the  time  that  would  be  required  to  bring  ' 

the  unit  to  steady  state  conditions,  the  time  constant  for  the  system  had 

to  be  estimated.  ,! 


The  thermal  losses  at  the  average  operating  temperature  of  1250°F  of  the 
thermal  energy  storage  unit  were  calculated  to  be  about  I80  watts.  These 

losses  are  slightly  higher  than  those  shown  in  Figure  21.  Thermal  losses  | 

from  the  fill  valve,  which  was  to  be  left  on  the  unit  until  the  operation  of 
the  unit  had  been  verified,  and  the  heat  losses  at  the  hot  cylinder  end  had 
to  be  added  to  the  calculated  losses.  The  effective  time  constant  was 
estimated  from  the  calculated  thermal  capacity  of  7.076  x 10^  j/°F  and  the 
predicted  thermal  losses  at  the  operating  temperature,  i.e.,  I80  watt  at 
T = 1250°F,  to  be  t = ^6.25  x 10^  sec  or  12.85  hours. 

This  indicated  that  the  unit  would  require  about  two  days  to  reach  its 
steady  state  temperature  at  each  test  point.  The  test  results  are  shown  in 
Figure  55.  Because  of  the  long  time  constant,  the  plotted  values  are  not 
absolutely  accurate,  as  it  would  have  required  a considerably  longer  time 
than  was  available  to  obtain  true  steady  state  data.  The  process  was 
speeded  up  by  establishing  the  temperature  loss  correlation,  estimating 
the  losses  at  a given  temperature,  and  adjusting  the  input  power  of  the 
heaters  to  equal  these  losses  when  a given  temperature  was  reached. 
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During  these  initial  tests,  one  of  the  three  heater  sets  was  used  to  supply 
the  steady  state  power  while  the  other  two  sets  wore  employed  to  raise  the 
temperature  level.  When  the  power  of  the  two  sets  of  heaters  was  raised 
for  an  average  power  output  per  heater  of  l62  watts  (this  compares  with  a 
design  power  output  of  672  watts),  three  heaters  failed  almost  simulta- 
neously. Inspection  of  the  heaters  indicated  that  all  three  heaters  had 
failed  in  the  very  same  mode.  Arcing  between  the  two  wires  of  each  heater 
had  occurred  at  the  location  where  the  transition  wire  had  been  welded  to 
the  lead  wire,  as  shown  in  the  sketch  below. 


TfiANSITION  LEAD  (NICKEL) 


The  heater  vendor  was  contacted.  Upon  reviewing  the  heater  design,  the 
responsible  manager  conceded  that  a "human  error"  had  occurred  in  the 
manfacturing  process.  It  was  felt  that  the  remaining  three  heaters  would 
not  survive  the  planned  testing,  and  that  it  would  be  appropriate  to 
replace  all  heaters  immediately  with  new  heaters. 


All  six  heaters  were  pulled  out  of  the  well  after  removing  the  insulation. 
New  heaters  from  a different  vendor  were  ordered  and  shipped  from  the  east. 
Although  the  new  heaters  had  the  same  power  capacity  as  the  original 
heaters,  their  power  loading  was  considerably  higher  than  that  of  the 
original  heaters  because  each  heatei’  was  only  29  inches  long.  This 
compared  with  the  70  inch  length  of  the  original  heaters.  Each  heater  was 
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bent  into  a hair  pin  and  pushed  into  the  heater  well  of  the  thermal  storage 
unit.  After  replacing  the  insulation,  testing  of  the  thermal  energy 
storage  unit  was  continued. 


During  testing  of  the  thermal  energy  storage  unit,  two  of  the  six  heaters 
of  the  second  set  of  heaters  failed  by  shorting  internally  to  a dead  short, 
although  remaining  isolated  from  the  sheath.  The  vendor  of  the  heaters  was 
contacted  who  explained  the  failure  by  an  accumulation  of  moisture  at  the 
cold  end  of  the  heater  (i.e.,  the  non-heated  extension  of  the  lead  wires) 
where  the  short  between  the  two  lead  wires  occurred.  He  asked  for  the 
return  of  all  six  heaters  for  repair  without  charge.  It  was  agreed  that 
the  heaters  would  be  returned  after  completion  of  testing  with  the  remain- 
ing four  heaters.  Since  each  heater  had  a capacity  of  over  1700  watts, 
ample  capacity  remained  for  continuing  testing.  The  steady  state  heat  loss 
data  were  found  not  to  vary  substantially  from  the  original  test  data; 
indicating  that  the  insulation  was  the  same  as  before.  The  heat  loss  could 
be  expressed  as  a function  of  temperature  by 

= A(AT)" 

where 

= power  loss,  watt 
A = 3.318  X 10“^ 

n = 1.5^16 

£.T  = temperature  difference  T - T , 

o room’ 

4. 2. 1.2  Thermal  Capacity 


With  the  correlation  of  the  thermal  losses  against  the  operating  tempera- 
ture known,  the  effective  thermal  capacity  of  the  thermal  energy  storage 
unit  could  be  determined  by  the  relation 


(p.  - A X (AT)"  - P . i 

V in  out  / 


dt 


dT 
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where 


= Effective  sensible  heat  of  the  entire  unit 


P. 

in 


out 


= Power  input  supplied  by  the  electric  heaters 

= Power  extracted  through  the  hot  cylinder  by  cooling  gas 
flow 


A,n  = Constants  determined  from  steady  state  energy  losses 
from  the  thermal  energy  storage  unit 
= Temperature  difference  between  thermal  energy  storage 
unit  temperature  and  ambient  temperature 


This  relation  was  computerized  to  permit  calculation  of  the  effective 
thermal  capacity,  thermal  losses,  and  power  output.  The  required  input 
data  were  the  electrical  power  input  and  operating  temperature  of  the  unit 
at  given  time  intervals. 


Below  and  above  the  eutectic  temperature  of  the  thermal  energy  storage 
material  the  effective  thermal  capacity  was  expected  to  be  relatively 
constant  showing  only  a slight  increase  with  increasing  temperature.  At 
the  eutectic  temperature,  the  effective  heat  capacity  should,  in  the  ideal 
case,  approach  infinity,  as  thermal  energy  is  either  stored  or  absorbed 
without  increase  in  temperature.  The  computer  program  was  expanded  to  also 
permit  calculation  of  the  energy  extracted  from  or  stored  by  the  thermal 
energy  storage  material  and  the  structural  material  of  the  thermal  energy 
storage  unit,  and  the  total  thermal  energy  extracted  as  a function  of  the 
temperature  or  the  thermal  energy  storage  unit. 

4.2.2  TEST  DATA  EVALUATION 


4.2.2. 1 Sensible  Heat  Capacity 

After  the  initial  tests,  which  provided  the  correlation  between  the 
temperature  of  the  thermal  energy  storage  unit  and  the  thermal  losses  from 
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the  unit,  the  unit  was  heated  to  a temperature  that  was  assumed  to  be  below 
the  eutectic  temperature  of  the  thermal  energy  storage  material.  The 
maximum  measured  heater  well  temperature  was  1303°F  with  a power  input  of 
214  watts.  Thus,  no  storage  of  thermal  energy  in  form  of  latent  heat  of 
fusion  should  have  occurred  and  all  energy  that  would  be  released  due  to 
thermal  losses  from  the  unit  through  the  insulation  would  be  the  total 
sensible  heat  of  all  materials  of  the  unit.  By  determining  the  temperature 

j of  the  unit  as  a function  of  time  and  applying  the  relation  shown  on  page  42 

it  should  be  possible  to  verify  the  calculated  thermal  capacity  of  the 
unit,  which  must  be  known  for  evaluating  the  thermal  capacity  of  the 
thermal  energy  storage  material. 

The  measured  correlations  between  thermal  energy  storage  temperature  and 
time  are  shown  in  Figures  66  and  67.  The  plotted  test  results  immediately 
indicated  that  the  temperatures  do  not  correlate  in  the  expected  manner. 

‘ Due  to  the  higher  thermal  losses  at  the  higher  temperatures,  a steeper  drop 

in  temperature  at  the  beginning  of  the  cooling  cycle  would  have  been 
predicted  if  the  unit  had  a relatively  constant  thermal  capacity.  On  the 
other  hand  the  data  indicated  a much  larger  thermal  capacity  at  the  higher 
temperature.  When  the  data  were  reduced,  the  thermal  capacity  was  found  to 
increase  with  temperature  as  shown  in  Figure  68. 

In  Figure  69,  the  temperatures  of  the  thermal  energy  storage  unit  are  shown 
as  they  were  measured  during  the  heating  of  the  unit.  The  appearance  of 
the  heating  curve  is  very  similar  to  the  appearance  of  the  cooling  curves, 
indicating  that  the  effective  thermal  capacity  of  the  unit  is  the  same 
whether  measured  during  cooling  or  during  heating. 

The  test  data  made  it  appear  that  the  effective  thermal  capacity  of  the 
thermal  energy  storage  unit  increased  with  increasing  temperature.  At  the 
average  operating  temperature  of  T = 650°F,  the  effective  thermal  capacity 
of  the  unit  appeared  to  he  7.22  x 10^  joules/°F.  This  compared  with  the 
predicted  effective  thermal  capacity  of  7.1  x 10^  joules/°F,  a variation  of 
1.7  percent.  The  abnormal  behavior  was  evident  from  the  cooling  curves  as 

( 

j 

i 
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well  as  from  the  heating  curves.  Despite  the  fact  that  no  melting  should 
have  occurred,  the  slope  of  the  curves  changed  their  trend  only  below 
950°F,  which  lies  considerably  below  the  reported  melting  point  of  the 
thermal  energy  storage  salt.  It  was  thought  that  the  variable  effective 
heat  capacity  of  the  thermal  energy  storage  unit  could  possibly  be 
explained  by  the  varying  amount  of  the  working  fluid  of  the  heat  pipe, 
sodium,  in  the  vapor  state. 

With  increasing  temperature,  the  amount  of  sodium  in  the  vapor  state 
increases.  The  effective  total  specific  heat  of  the  sodium  vapor  due  to 
evaporation  and/or  condensation  in  the  vapor  volume  is  given  by  the  change 
of  latent  heat  capacity  with  temperature,  i.e., 


S^eff  Na  = ^ 


where 


AHxp  xVxM 
*^v 


latent  heat  of  evaporation  content  of  sodium  vapor  in 

the  vapor  volume,  joules 

latent  heat  of  evaporation,  joules/gram 

vapor  pressure,  dynes/cm^ 

3 

volume,  cm 


= molecular  weight,  grams/mole 
o 

= temperature,  K 
= gas  constant,  ergs/°K-mole 


The  vapor  pressure  of  sodium  can  be  expressed  as  a function  of  temperature 
by  the  relation 


Py  = A exp  (-E/T) 


I 


where 


10  2 
3.363  X 10  dynes/cm 

11973. 'I^K 


The  latent  heat  of  evaporization  can  be  found  by  employing  Clapeyron's 
relation 

AH  = T(v^  - v^)  dp^/dT 


where 


= vapor  volume 
v^  = fluid  volume 


Since 


V » V 
1 2 


the  relation  can  be  reduced  to 


AH  = T X V,  X dp  /dT 
1 ^v 


By  differentiating  relation  3 and  employing  the  relation 


p X V = RT/M 


the  latent  heat  of  evaporization  can  be  found  to  be 


AH  = RE/M 


The  effective  thermal  capacity  of  the  sodium  can  then  be  expressed  by 


'p  'eff  Na 


E X (H/T  - 1)  x V X A X exp  (-F./T) 


r 


The  total  volume  V,  of  the  vapor  in  the  heat  pipe  of  the  unit  was  calculated 
to  be  1243.9  cm^. 

When  the  above  relation  was  evaluated,  as  shown  in  the  following  tabula- 
tion, the  change  of  state  of  the  sodium  was  found  to  have  only  a minor 
effect  on  the  apparent  thermal  capacity  of  the  entire  system.  It  could  not 
explain  the  observed  change  in  thermal  capacity  of  the  system. 

T,°K  300  500  750  1000 

CpW,j/°K  1 X 10'^°  1.8  X 10"*^  0.16  3.48 

4. 2. 2. 2 Latent  Heat  Capacity 

After  the  completion  of  these  initial  tests,  the  temperature  of  the  thermal 
energy  storage  unit  was  increased  above  the  melting  temperature  of  the 
thermal  energy  storage  salt.  The  electric  power  to  the  heater  was  then 
turned  off  and  the  unit  was  permitted  to  cool  solely  by  the  losses  through 
the  insulation.  The  results  of  that  test  are  shown  in  Figures  70  and  71. 
It  can  be  seen  that  the  solidification  of  the  thermal  energy  storage  salt 
started  when  the  temperature  of  the  thermocouple  on  the  outside  of  the 
thermal  energy  storage  unit  reached  1280°F. 

When  the  data  were  reduced  to  determine  the  amount  of  thermal  energy 
extracted  from  the  unit  as  function  of  time  and  temperature,  the  distribu- 
tion of  the  total  energy  between  sensible  heat  and  phase  change  energy  as 
shown  in  Figures  12  and  73  wore  found.  The  results  seemed  to  indicate  that 
the  total  phase  change  energy,  which  had  be»n  expected  to  consist  primarily 
of  the  fusion  energy  of  the  thermal  energy  storage  material,  was  not 
released  within  a -25*^F  temperature  range. 

The  thermal  energy  storage  unit  design  called  for  a thermal  energy  storage 
of  4.32  X 10^  joules  in  the  form  of  fusion  energy.  The  amount  of  thermal 
energy  storage  material  required  was  based  on  a fusion  energy  of 
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Thermal  Energy  Storage  Demonstration  Unit  During  Cool-down 


S6729 


TEMPERATURE  T.  ®F 

Figure  72.  Extracted  Energy  (Cooling  by  Thermal  Losses  Only) 


788.4  joules/gram  (339.27  B/lb  or  0.219  kW-hr/kg) . This  compares  with  the 
fusion  energy  of  944  joules/gram  based  on  ideality  of  the  ternary 
eutectic  mixture  of  MgF2~LiF-KF.  Based  on  ideality,  the  total  fusion 
energy  available  in  the  thermal  energy  storage  unit  would  have  been 
5.1?  X 10^  joules. 

From  the  results  plotted  in  Figures  72  and  73  it  appears  that  only 
2.45  X 10^  joules  of  phase  change  energy  was  released  over  a temperature 
range  of  50°F.  This  is  about  56.7  percent  of  the  design  energy.  The  total 
amount  of  the  design  energy  was  extracted  over  a temperature  range  of  1285° 
to  1055°F,  i.e.,  230°F.  The  total  amount  of  phase  change  energy  appeared 
to  be  5.07  X 10^  joules.  This  is  almost  equal  to  the  amount  of  energy  that 
would  have  been  calculated  to  be  available  as  fusion  energy  based  on 
ideality. 

The  same  energy  content  in  the  TES  salt  as  a function  of  temperature  is 
exhibited  in  the  heating  curve  as  shown  in  Figure  74.  This  verifies  that 
the  energy  content  is  independent  of  either  the  cooling  or  heating  process 
of  the  salt. 

The  effective  thermal  capacity  of  the  thermal  energy  storage  unit  as  a 
function  of  time  is  shown  in  Figure  75.  The  curve  shows  two  peaks  in  the 
effective  thermal  capacity  of  TES  unit,  one  at  about  1290°F  and  the  other 
at  1165°F.  At  these  two  temperatures  thermal  energy  is  released  with  the 
smallest  drop  in  temperature.  It  is  not  known  what  phase  change  coincides 
with  the  second  peak  at  1165°F. 

4.2.2. 3 Effect  of  Extraction  Rate 

The  first  evaluation  of  the  thermal  energy  storage  capacity  was  from  test 
data  obtained  during  a test  in  which  the  thermal  energy  storage  unit  was 
cooled  only  by  its  losses,  which  were  a function  of  the  temperature  of  the 
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Figure  75.  Apparent  Ther.Tial  Capacity  of 
Thermal  Energy  Storage  Unit 
During  Charging 


A 


unit.  The  initial  loss  was  213  watts  decreasing  with  decreasing  tempera- 
ture, but  remaining  almost  constant  during  the  time  at  which  the  majority 
of  the  fusion  energy  was  released. 

Test  results  are  shown  in  figures  76  and  77  of  a test  at  which  the  initiai 
total  thermal  power  extraction  from  the  unit  was  555.6  watts.  The  power 
was  extracted  from  the  unit  by  losses  througii  the  insulation  and  by  cooling 
the  hot  cylinder  by  an  air  flow  of  W = 0.76  gram/sec.  At  the  eutectic  point 
the  power  extraction  rate  was  475  watts.  In  a third  test  whose  test  data 
are  presented  in  Figures  78  and  79,  the  gas  flow  rate  was  W = 3-7  gram/sec. 
At  this  flow  rate  the  initial  power  extraction  rate  was  1799  watts,  while 
at  the  eutectic  temperature  the  extraction  rate  was  1067  watts,  which  was 
almost  equal  to  the  design  extraction  rate  of  1200  watts,  200  watts  of 
losses  and  1000  watts  of  power  carried  away  by  the  cooling  gas  in  the  hot 
cylinder. 

When  the  amount  of  energy  extracted  from  the  thermal  energy  storag'-’ 
material  is  correlated  with  the  extraction  rate,  no  effect  of  the  extrac- 
tion rate  on  the  evaluated  amount  of  available  thermal  energy  can  be  found. 
The  variations  in  the  reduced  data  are  within  the  accuracy  with  which  the 
temperatures  and  the  times  could  be  measured. 

Thus,  the  test  data  indicated  that  the  behavior  of  the  thermal  energy 
storage  material  selected  for  this  application  was  different  from  that 
anticipated.  While  lithium  fluoride  liberated  the  entire  fusion  energy  at 
a specific  melting  point,  the  ternary  eutectic  salt  made  up  of  MgF^,  LiF, 
and  KF,  liberates  phase  change  energy  over  a wide  range  of  temperature. 
Since  this  temperature  range  was  not  known  at  the  time  of  the  design  of  the 
thermal  energy  storage  unit,  it  was  not  included  in  the  calculation  of  the 
temperature  swing  over  v;hich  the  hot  cylinder  would  operate.  Only  the 
increasing  temperature  drop  associated  with  the  build-up  of  solidified 
thermal  energ.y  storage  material  had  been  considered.  The  overall  effect  of 
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Figure  77.  Extracted  Energy  (Cooling  by  Therynal  Los.'so.'?  and  by 
Simulated  V/M  Cooler  Operation  with  Flowing  Gas) 
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the  unpredicted  behavior  of  the  chosen  TES  material  is  either  the  accep- 
tance of  a larger  temperature  swing  over  which  the  hot  cylinder  of  the 
Vuilleumier  cooler  has  to  operate  or  a smaller  energy  storage  density. 
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4. 2. 2. 4 Temperature  Distribution 

The  design  of  the  thermal  energy  storage  unit  was  primarily  aimed  at 
extracting  energy  from  the  thermal  energy  storage  material  and 
transferring  it  to  the  hot  cylinder  with  the  smallest  possible  temperature 
drop.  This  would  permit  the  use  of  a thermal  energy  storage  material  that 
had  a melting  temperature  only  slightly  higher  than  the  operating  tempera- 
ture of  the  hot  cylinder  of  the  Vuilleumier  cooler.  The  best  approach  for 
achieving  the  desired  operating  characteristic  of  the  thermal  energy 
storage  unit  was  the  use  of  the  heat  pipe  phenomenon  by  which  energy  could 
be  transferred  over  a relatively  large  distance  with  a relatively  small 
temperature  drop. 

When  the  thermal  energy  storage  unit  was  tested  with  energy  extraction  at 
the  hot  cylinder,  large  temperature  gradients  were  measured  in  the  thermal 
energy  storage  unit.  In  Figures  80  and  8l  the  temperatures  which  were 
measured  and  recorded  on  a strip  chart  and  whose  locations  are  presented  in 
Figure  60  are  plotted  against  time.  The  gas  flow  rate  through  the  hot 
cylinder  was  set  to  extract  a total  power  of  1067  watts.  Since  the 
temperature  differences  did  not  exist  during  charging,  as  can  be  seen  from 
Figure  82,  the  observed  temperature  gradients  could  not  be  attributed  to 
faulty  thermocouple  readings. 

All  evidence  indicated  that  the  heat  pipe  functioned  only  partially.  Some 
heat  pipe  effect  was  deducted  from  the  fact,  that  the  power  transferred  at 
the  hot  cylinder  could  not  be  conducted  through  the  structural  components 
of  the  the  thermal  energy  storage  unit  with  the  indicated  temperature 
gradients  along  the  unit. 


113 


Three  possible  reasons  for  the  malfunction  of  the  heat  pipe  could  be 
considered;  (1)  insufficient  amount  of  working  fluid  (sodium),  (2)  failure 
of  the  wick,  and  (3)  incondensable  gas  in  the  heat  pipe.  All  three 
possible  modes  of  heat  pipe  deficiency  were  investigated.  (See  Appendix  B 
for  wick  study.)  The  failure  was  finally  traced  to  a gas  leak  in  the  fill 
valve. 

The  heat  pipe  was  reprocessed  after  the  fill  valve  had  been  replaced  by  a 
new  valve.  When  the  stored  energy  was  extracted  from  the  unit  by  thermal 
losses  *nd  by  cooling  of  the  hot  cylinder,  the  thermal  energy  storage  unit 
exhibited  the  temperature  distribution  that  had  been  predicted  for  the 
design  (Figure  83).  All  temperatures  along  the  outside  of  the  unit 
remained  equal,  while  the  temperature  measured  inside  the  thermal  energy 
storage  unit  (T.C.  Nos.  1 and  2)  was  higher.  The  difference  between  the 
outside  and  inside  temperatures  was  increasing  with  time  as  had  been 
expected,  when  considering  the  heat  transfer  across  the  solidified  TES 
material.  The  average  power  extraction  rate  was  about  883  watts  (fig- 
ure 84)  and  the  temperature  differential  between  the  outside  temperatures 
and  the  inside  temperature  was  42°F.  A temperature  differential  of  50°F  at 
an  extraction  rate  of  1200°  watt  had  been  predicted  by  the  heat  transfer 
analysis  that  was  presented  in  Subsection  2.2.1. 

The  measured  energy  extracted  as  a function  of  temperature  is  shown  in 
Figure  85. 


116 


SOO  600  700 


aoo  900  1000 

temperature  t,  “f 


Figure  85.  Extractpr)  Enorfr.y  (Coolinr,  by  Therrml  Lorir.or.  and  by 
Simulalod  V/M  Cooler  Operation  with  Flowing  Gan) 
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SECTION  V 


CONCEPTUAL  FLIGHT  DEVICE  PARAMETER  DEVELOPMENT 


5. 1 BASIC  REQUIRMENTS 

For  the  evaluation  of  the  thermal  energy  storage  concept  as  applied  to  a 
high  capacity  Vuilleumier  refrigerator,  requirements  for  the  refrigerator 
as  shown  in  Table  6 were  assumed. 

TABLE  6.  REFRIGERATOR  REQUIREMENTS 


Orbit : 

12  hr 

Eclipse  Time: 

1 +0.1  hr 

Operating  Temperature: 

1250°  +25' 

Total  Power  Requirement: 

1000W 

1500W 

2000W 

No.  of  Hot  Cylinders: 

1 or  2 

With  the  Vuilleumier  cooler  energy  requirement  given,  it  was  further 
assumed  that  thermal  losses  would  amount  to  about  10  percent  of  the  cooler 
power,  and  that  the  thermal  energy  storage  material  is  a ternary  eutectic 
salt  made  up  of  64  LiF  - 30  MgF^  - 6 KF,  having  an  eutectic  temperature  of 
1310°F  and  a heat  of  fusion  of  788.4  joules/gram  (339.3  B/lb). 

The  main  restrictions  placed  on  the  design  of  a thermal  energy  storage  unit 
were: 

I 

I 

a.  The  unit  must  be  constructed  from  commercially  available  Inconel 

600  pipes.  j 

i 

1 19 




b.  The  unit  must  be  tested  in  a 1-g  environment. 

c.  The  unit  will  be  mrited  to  the  hot  cylinder  of  the  High  Capacity 
refrigerator. 

The  stipulation  that  the  design  should  incorporate  commercially  available 
Inconel  600  pipes  was  based  on  the  experience  gained  in  trying  to  procure 
specially  made  tubing  from  sheet  metal  stock.  An  intensive  procurement 
inquiry  produced  only  two  vendors  who  were  willing  to  undertake  the 
fabrication  of  some  but  not  all  tubes  to  a given  design.  Despite  quoting  a 
price  three  times  higher  than  that  quoted  by  the  first  vendor,  the  second 
vendor  was  not  even  willing  to  guarantee  the  loosest  tolerances  required 
for  final  welding  of  the  thermal  energy  storage  unit.  Thus,  the  available 
sources  for  tube  fabrication  from  sheet  metal  stock  proved  to  be  too  few 
and  the  fabricatior  cost  appeared  to  be  relatively  high. 


The  nonfunctional  requirement  imposed  on  the  design  that  the  unit  be  able 
to  operate  in  a 1-g  environment,  has  to  be  accepted  as  long  as  the  therril 
energy  storage  unit  has  to  be  tested  and  its  operation  verified  prior  to 
launching  and  operation  in  space.  This  requirement  imposes  a dimensional 
constraint  on  the  unit  which  is  set  by  the  height  to  which  the  working 
fluid  of  the  heat  pipe  can  be  raised  in  a 1-g  environment  by  the  suction 
pressure  due  to  capillary  forces  in  the  wick  structure  of  the  heat  pipe. 


5.2  CONCEPTUAL  DESIGN  PARAMETERS 


5.2.1  BASIC  DESIGN 


The  basic  design  that  evolved  from  the  above  requirements  and  constraints 
is  shown  in  Figure  86.  The  length,  L,  of  the  thermal  energy  storage  unit  is 
determined  by  the  power  requirement  of  the  refrigerator  which  is  shown  in 
Figure  87. 
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Figure  87.  Corrolntion  Between  Power  Requirement 
and  TF,S  Unit  i.ength 
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5.2.2  TEMPERATURE  DISTRIBUTIONS 

The  temperature  distribution  in  the  thermal  energy  storage  unit  is  sho'.m  in 
Figure  88.  This  temperature  distribution  is  not  affected  by  the  length  of 
the  unit,  because  the  temperature  drop  associated  with  the  vapor  flow 
pressure  drop  is  very  much  less  than  I'^F.  It  remains  negligible  over  the 
entire  range  of  design  lengths  considered  in  this  study. 

The  temperature  drop  across  the  thermal  energy  storage  material  as  a 
function  of  the  discharge  time  is  shown  in  Figure  89.  The  effect  of  the 
temperature  drop  across  the  thermal  energy  storage  material  on  the  operat- 
ing temperature  of  the  Vuilleumier  refrigerator  is  shown  in  Figure  90.  At 
the  beginning  of  the  discharge  cycle,  the  operating  temperature  is  1279°F, 
while  at  the  end  of  the  cycle  the  operating  temperature  will  be  either 
1227°F  or  1219°F  for  a 1-hour  orbit,  respectively.  The  operating  tempera- 
ture is  therefore  29°F  above  the  design  operating  temperature  of  at 
the  beginning  of  the  eclipse  and  23°F  below  the  design  temperature. 

The  operating  temperature  of  the  Vuilleumier  refrigerator  over  the  entire 
12-hour  cycle,  which  includes  the  charging  and  discharging  of  the  thermal 
energy  storage  unit,  is  shown  in  Figure  91.  It  would  be  interesting  to 
compare  this  temperature  variation  and  the  cycle  frequency  of  a hot 
cylinder  operated  with  a thermal  energy  storage  unit  with  the  temperature 
history  of  a hot  cylinder  which  is  heated  directly  with  a temperature 
controlled  electrical  heater. 

5.2.3  WEIGHT 

The  total  weight  of  the  thermal  energy  storage  unit  including  all  necessary 
components  of  the  unit,  structural,  wicking,  working  material,  and  insula- 
tion, i.s  presented  in  Figure  92.  For  comparison,  the  probable  weight  of  a 
nickel-hydrogen  battery  that  can  satisfy  the  same  power  requirements  as  the 
th*'rmal  en'^r'gy  storag'’  unit  are  presented.  The  specific  weight  of  the 
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Figure  83.  Temperature  Distribution  in  the  Thermal  Energy  Storage  Unit  and  the  Hot 
Cylinder  of  the  Long  Life,  High  Capacity  Vuilleumier  Refrigerator  for 
Space  Application 
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thermal  energy  storage  unit  is  shown  in  Figure  93-  With  increasing  power 
requirements,  the  specific  weight  approaches  asymtotically  3^  watt-hr/lb 
for  a single  hot  cylinder  Vuilleumier  refrigerator,  and  slightly  more  than 
32  watt-hr/lb  for  a dual  hot  cylinder  refrigerator.  The  specific  weights 
for  batteries  are  being  presented  only  for  initial  evaluation. 

When  comparing  the  weights  of  the  thermal  energy  storage  material 
(Figure  94)  with  the  weights  of  the  container  (Figure  95),  it  can  be  seen 
that  they  are  equal.  While  the  thermal  energy  storage  material  weight  is 
determined  by  the  enei-gy  storage  requirement,  the  container  weight  could 
conceivably  be  reduced  if  and  when,  by  experience  and  long  duration 
testing,  the  minimum  material  thickness  for  the  container  material  has  been 
established, 

5.2.4  STRUCTURE 

For  the  outer  container,  the  wall  thickness  is  determined  by  the  collapsing 
criteria  of  the  cylinder.  The  unit  has  to  withstand  a pressure  difference 
of  one  atmosphere  when  the  thermal  energy  storage  \mit  is  cold  and  a 
slightly  lower  pressure  differential  when  the  unit  is  operating  at  its 
design  temperature.  Correspondingly,  the  materials  properties  of  Inconel 
600  are  changing  with  temperature  as  shown  in  Figure  96.  The  correlation 
between  wall  thickness  of  the  outer  cylinder  and  the  collapsing  pressure  is 
shown  in  Figure  97.  If  a safety  factor  of  4 is  employed,  then  the  wall 
thickness  for  the  design  indicated  in  Figure  86  would  be  about  0.060  inch. 

The  collapsing  pressure  calculations  assumed  a long  cylinder  to  cover  the 
entire  range  of  design.  The  calculations  can  be  refined  to  incorporate  the 
effect  of  length  on  the  wall  thickness  requirements.  A relatively  short 
TES  unit  could  be  constructed  with  a slightly  thinner  wall  than  that 
indicated  by  Figure  97. 
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Figure  93.  Specific  Energy  of  n Thermal  Energy  Storage  Unit  for  a Vui  1 1 e i-.ier 
Cooler 
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Figure  96.  Dynamic  Modulus  of  Elasticity  E and  Yield  Strength 
of  Inconel  600  Alloy 
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Figure  97.  Collapsing  Pressure  of  Thin- 
Walled  Vessel 
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5.2.5  CONCEPTUAL  FLIGHT  DEVICE 

For  a mission  using  a long  life,  high  capacity  Vuilleumier  refrigerator 
with  two  opposing  hot  cylinders,  the  total  thermal  power  requirement  for 
the  refrigerator  is  1150  watts.  This  power  has  to  be  supplied  to  the  two 
opposing  hot  cylinders,  i.e.,  575  watts  to  each  cylinder.  The  low  earth 
orbit  would  demand  only  18  minutes  of  thermal  energy  to  be  supplied  by  the 
TES  units.  Using  the  design  underlying  the  TES  demonstration  unit  which 
called  for  an  energy  requirement  of  1000  watts  for  60  minutes,  the  length 
of  the  TES  units  for  the  long  life,  high  capacity  Vuilleumier  refrigerator 
(Hi-Cap  V/M)  hot  cylinders  can  be  estimated  to  be 

T - ^^50  ^28 

■ 2 X 1000  60  TESDU  " TESDU 

The  hot  cylinder  physical  layout  with  the  TES  unit  is  shown  in  Figure  98 
and  Figure  99  shows  the  Hi-Cap  V/M  cooler  with  the  TES  units  mounted  on 
both  hot  cylinders.  For  comparison.  Figure  100  shows  the  cooler  without 
the  TES  units. 

The  available  pipe  sizes  for  the  construction  of  the  TES  unit  are  boxed  in 
black  in  Figure  32.  Only  Schedule  40  pipes  are  available  for  larger  pipe 
sizes.  The  nominal  5-inch  pipe  has  a wall  thickness  of  0.258  inch,  which 
is  about  three  times  thicker  than  needed.  Thus,  an  unraachined  pipe  would 
weigh  about  10  pounds,  but  approximately  two-thirds  (6.5  pounds)  would  have 
to  removed.  The  cost  of  the  pipe  is  approximately  $l6/pound  indicating  a 
waste  of  about  $100. 


Figure  99.  Long  Life,  High  Capacity  Vuilleumier  Refrigerator  with 
Thermal  Energy  Storage  Units 


Figure  100.  Lon‘t  I.ife,  High  Capacity  Vuilleumier  Refrigerator  with 
Battery  Fnerg.ir.eh  Kieet.ric  Heaters  (Present  Hi-Cap 

Configuration) 


SECTION  VI 


CONCLUSIONS 


The  Thermal  Energy  Storage  Demonstration  Unit  program  had  as  its  objective 
proof  of  feasibility  of  storing  energy  in  form  of  latent  heat  for  powering 
a Vuilleumier  cooler  during  the  time  a satellite  is  in  the  eclipse  of  the 
sun.  The  major  technical  areas  that  needed  technical  confirmation  were  the 
storing  of  thermal  energy,  the  extraction  of  the  energy  over  a limited 
temperature  range,  and  the  transfer  of  the  thermal  energy  to  the  hot 
cylinder  of  the  Vuilleumier  cooler. 

The  design  of  the  thermal  energy  storage  unit  confirmed  the  concept  of 
storing  energy  from  an  electric  heater  in  form  of  the  latent  heat  of  a 
thermal  energy  storage  material.  The  stored  thermal  energy  was  extracted 
over  a large  area  with  a minimum  of  temperature  drop  and  transferred  to  the 
hot  cylinder  by  heat  pipe  action. 

The  testing  of  the  thermal  energy  storage  demonstration  unit  reveal<’d 
rather  unexpectedly  that  the  phase  change  energy  of  the  ternary  eutectic 
salt,  LiF-MgF2“KF,  was  not  released  or  absorbed  entirely  at  the  melting 
point  of  the  salt,  but  over  a relatively  wide  temperature  range.  It 
appeared  that  only  about  50  percent  of  the  "ideal"  fusion  energy  is 
absorbed  during  hf'ating  or  r<^leaned  during  cooling  at  th<^  eutectic  tempera- 
ture of  about  1310°F.  Phase  tranr.formation  with  energy  absorption  upon 
heating  apparently  takes  place  at  t''mperatuiros  below  the  eutec'ic  tem.pera- 
ture.  This  pro'diic'^s  ‘he  appearan''<'  of  a specific  heat  which  is  m.uch  m.ore 
temperature  lepi^n  lent  •'han  kinetic  ‘.heorv  would  predict. 

The  thermal  enerv  • »' i.u-  b-m  Vi  ■ * • u'.  u-:  ‘ pr>  gram  demonstrated  a TF.S 
unit  that  store]  p.,  , ^ enet'gy,  released  the  stored 


1 I' 


energy  at  the  eutectic  temperature  and  below,  and  could  transfer  the  energy 
by  heat  pipe  action  to  the  hot  cylinder  of  a Vuilleumier  cooler.  The  large 
amount  of  ther'mal  energy  released  at  the  eutectic  temperature  and  during 
the  solidification  of  the  thermal  energy  storage  salt  was  transferred  over 
a temperature  range  of  about  This  agreed  well  with  the  design 
prediction,  confirming  the  capability  of  designing  TES  units  to  specific 
requirements  when  the  thermodynamic  properties  of  the  TES  material  are 
known . 

The  evaluation  of  weight  contributions  of  the  various  TES  unit  components 
pointed  to  the  most  fruitful  areas  of  further  investigation  of  all  aspects 
of  thermal  energy  storage,  if  this  concept  should  prove  to  be  a valid 
alternative  to  storing  energy  in  an  electrical  battery.  A considerable 
reduction  of  the  wall  thicknesses  of  the  container  could  be  possible,  if 
material  compatibility  for  long  duration  operation  could  be  proven  to  exist 
between  the  storage  salt,  the  heat  pipe  working  fluid,  and  the  container 
material . 
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EFFECT  OF  NONCONDENSABLE  GAS  ON  HEAT  PIPE  OPERATION 


When  large  temperature  differences  along  the  heat  pipe  section  of  the 
thermal  energy  storage  unit  were  observed,  an  attempt  was  made  to  explain 
this  apparent  malfunction  of  the  heat  pipe  by  assuming  the  possibility  of 
noncondensable  gas  in  the  heat  pipe  section  of  the  thermal  energy  storage 
unit.  In  Figure  A-1,  the  temperature  difference  between  thermocouple  No.  5 
and  thermocouple  No.  10  were  plotted.  The  temperature  difference  could  be 
explained  by  a vapor  pressure  difference  of  the  working  fluid  which  has 
been  calculated  for  the  operating  temperature  and  which  has  also  been 
plotted  in  Figure  A-1.  The  calculated  pressure  di f ferentials  which  are 
based  on  the  observed  temperature  dif ferentials,  ar«  considerably 
higher  than  those  which  were  calculated  for  the  vapor  flow,  as  can  be  seen 
from  the  pressure  differentials  shown  in  Figures  l6,  17  and  18  of  Subsec- 
tion 2. 2. 2. 3.  Only  the  pre.sence  of  noncondensable  gas  in  the  heat  pipe 
could  explain  such  largo  vapor  pressure  differentials.  If  the  presence  of 
noncondensable  gas  at  a pressure  of  1.6?  mm  Hg  at  room  temperature  is 
as.sumed,  then  temperature  differentials  as  shown  in  Figure  A-1  could  be 
expected.  The  calculated  temperature  differential  curve  has  apparently  a 
different  tr(!nd  than  the  curve  which  was  constructed  from  the  measured 
temperature  di  fferentials . Nevertheless,  some  similarity  in  b>’havii''' 
seem;;  to  exist. 


It  was  interest.ing  to  show  that,  if  noncondcnsable  gas  should  have  an 
effect  on  the  heat  pipe  operation,  a heat  pipe  operating  with  potassium  as 
the  working  fluid  would  he  less  affected  than  a heat  pipe  opei’ating  with 
sodium  at  the  same  op'^rating  temperature.  This  can  be  seen  in  Figure  A-2. 
The  reason  for  ttiis  diffeence  between  the  two  working  fluids  can  be  foun 1 
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from  the  graph  showing  the  vapor  pressure  of  sodium  and  potassium 
(Figure  A-3).  For  the  thermal  energy  storage  unit  operating  at  around 
1300  F,  potassium  would  be  quite  acceptable  as  heat  pipe  working  fluid,  as 
the  vapor  pressure  at  that  temperature  is  only  slightly  above  1/2  atm.  The 
unit  has  to  be  designed  for  a pressure  differential  of  1 atm  to  satisfy  the 
condition  in  the  laboratory. 
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TRANSIENT  WICK  EFFECTS 


When  initially  it  was  not  possible  to  verify  the  presense  of  noncondensable 
gas  in  the  heat  pipe  section  of  the  thermal  energy  storage  unit,  and  it 
appeared  that  the  heat  pipe  section  had  been  fully  evacuated,  the  capacity 
and  performance  of  the  wick  structure  was  investigated  in  an  informative 
experiment  to  verify  the  wick  design. 

The  wick  of  the  heat  pipe  was  constructed  with  a 70  Mesh  bolting  cloth 
which  has  a wire  diameter  of  0.00375  inch.  For  a fluid  the  wicking  height 
is 


h = A p^/p  g 

where: 

h = height,  cm 

2 

p = specific  gravity  of  water,  g/cm 

2 

g = gravitational  constant,  981  cm/sec 

A p^  = capillary  pressure,  dynes/cm 

A p = 2c  cosc/r 

“^c  ’ c 


where: 

c = surface  tension  of  water  against  air,  dynes/cm 
C*  = contact  angle,  degrees 
r^  = effective  pore  radius,  cm 

r = 0.707  ( 1/M  - d) 

c 
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where: 


M = Mesh,  cm 
d = wire  diameter,  cm 

For  distilled  water  at  room  temperature,  the  surface  tension  is 
c = 72.75  dynes/cm, 

and  the  wicking  height  is  found  to  be: 
h : 3-086  inch. 

Figure  B-1  shows  a fixture  that  was  built  to  simulate  the  wick  structure  in 
the  heat  pipe  section  of  the  thermal  energy  storage  unit.  In  step  one,  the 
dry  structure  was  weighted.  In  step  two,  the  wick  structure  was  soaked 
with  water  and  the  structure  was  weighted  again.  In  step  three,  the 
structure  was  held  horizontally  in  water,  after  which  it  was  weighted 
again.  In  the  fourth  step,  the  wick  was  permitted  to  soak  up  water  after 
most  of  the  water  had  been  shaken  out  of  the  wick.  The  structure  was 
weighted  again.  The  fully  filled  '1.9  inch  high  screen  appeared  to  hold 
7.0  grams  of  water.  When  standing  vertically,  it  retained  4.5  grams  of 
water.  This  indicated  that  the  height  of  the  water  in  the  wick  was 

h = h X W/W 
o o 

h = 3.15  inch 

This  compared  with  the  predicted  value  of  h = 3-086  inch.  After  the  fourth 
step,  the  water  weight  was  W = 3-6  grams,  whicii  corresponded  to  a wieking 
height  of  h = 2.52  inch,  or  about  80  percent  of  the  maximum  wieking 
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height.  The  data  of  Ref.  3 are  difficult  to  correlate  and  de  not.  [)crtain  to 
sodium.  Nevertheless  it  is  indicated  that  time  is  a stron-'  parameter  which 
determines  the  height  during  continuous  operation  to  which  the  wick  will 
remain  wet  in  a heat  pipe  whose  diameter  is  close  to  the  calculated  maximum 
wicking  height. 

The  small  informative  test  indicated  that  some  basic  data  should  be 
generated  in  support  of  the  design  of  full  scale  thermal  energy  storare 
units.  Wick  data  in  the  literature  and  of  other  investigators  are  very 
speculative  and  sketchy  and  can  therefore  not  be  used  for  designs  in  which 
the  limit  of  the  capability  of  wick  structures  is  being  approached.  The 
permeability  investigation  of  wicking  structures,  which  was  performed 
under  the  previous  program,  proved  the  need  for  such  basic  studies  in 
connection  with  the  design  of  a major  hardware  component. 

Tests  might  be  needed  to  ascertain  whether  the  heat  pipe  section  of  the 
thermal  energy  storage  unit  is  not  subject  to  the  fast  start-up  phenomenon 
which  was  observed  with  the  15-ft  long  heat  pipe.  At  an  extraction  rate  of 
1200  watts,  0.2772  g/sec  of  sodium  has  to  circulate.  This  is  equivalent  to 
997.92  g/hr  of  sodium.  The  entire  heat  pipe  section  contains  329  grams  of 
sodium,  indicating  that  the  sodium  has  to  be  wicked  back  three  times  an 
hour  during  operation.  This  rate  might  make  the  heat  pipe  operation 
marginal  for  a wick  height  which  is  more  than  75  percent  of  the  maximum 
wicking  height.  Some  basic  investigations  with  water  should  be  able  to 
answer  these  questions. 


3-  H.  R.  Knuz,  L.  S.  Langston,  B.  H.  Hilton,  S.  S.  Wyde,  and  G.  H. 
Nashick,  "Vapor-Ghami)er  Fin  .Studies  - Transport  Properties  and  Boiling 
Characteristics  of  Wicks,"  NASA  Contractor  Report  No.  NA.IA  CR-812, 
dune  1967. 
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